Technical  Report  EL-95-34 
December  1995 


Dbdif 

DJlll 


US  Army  Corps 
of  Engineers 

Waterways  Experiment 
Station 


Effect  of  Soil  Composition 
on  Complex  Dielectric  Properties 

by  John  O.  Curtis,  Charles  A.  Weiss,  Jr.,  Joel  B.  Everett 


Approved  For  Public  Release;  Distribution  Is  Unlimited 


19960229  089 


Prepared  for  U.S.  Army  Corps  of  Engineers  Research 
and  Development  Directorate 

DTIC  QUALITY  IMSPBOTED  8 


The  contents  of  this  report  are  not  to  be  used  for  advertising, 
publication,  or  promotional  purposes.  Citation  of  trade  names 
does  not  constitute  an  official  endorsement  or  approval  of  the  use 
of  such  commercial  products. 


PRINTED  ON  RECYCLED  PAPER 


Technical  Report  EL-95-34 
December  1995 


Effect  of  Soil  Composition 
on  Complex  Dielectric  Properties 

by  John  O.  Curtis,  Charles  A.  Weiss,  Jr.,  Joel  B.  Everett 

U.S.  Army  Corps  of  Engineers 
Waterways  Experiment  Station 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 


Final  report 

Approved  for  public  release;  distribution  is  unlimited 


Prepared  for  U.S.  Army  Corps  of  Engineers  Research 
and  Development  Directorate 
Washington,  DC  20314-1000 


Waterways  Experiment  Station  Cataloging-in-Publication  Data 

Curtis,  John  O.  o  •  /-u  i 

Effect  of  soil  composition  on  complex  dielectric  properties  /  by  John  O.  Curtis,  Charles 

A.  Weiss,  Jr.,  Joel  B.  Everett ;  prepared  for  U.S.  Army  Corps  of  Engineers  Research  and 
Development  Directorate. 

283  p.  :  ill. ;  28  cm.  --  (Technical  report ;  EL-95-34) 

Includes  bibliographic  references.  .  i 

1.  Dielectric  measurements.  2.  Soils  —  Composition.  3.  Soils  —  Testing.  I.  Weiss,  Char  es 
A.  II.  Everett,  Joel  B.  III.  United  States.  Army.  Corps  of  Engineers.  IV.  U.S.  .^y 
Engineer  Waterways  Experiment  Station.  V.  Environmental  Laboratory  (U.S.  Army  Engineer 
Waterways  Experiment  Station)  VI.  United  States.  Army.  Corps  of  Engineers.  Research  and 
Development  Directory.  VII.  Title.  VIII.  Series:  Technical  report  (U.S.  Army  Engineer 
Waterways  Experiment  Station) ;  EL-95-34. 

TA7  W34  no.EL-95-34 


Contents 


Preface  . . .  vi 

1—  Introduction  . 1 

Motivation  for  Study . 1 

Supporting  Studies  . ..2 

Fundamental  Relationships  . 3 

Microwave  Loss  Mechanisms  . 6 

Ionic  conductivity . 6 

Free  water  relaxation . 6 

Bound  water  relaxation  . 8 

Maxwell-Wagner  effect  . 8 

Study  Objectives  . 8 

2 —  Experimental  Technique  for  Measuring  Complex  Dielectric 

Constants . , . 9 

Governing  Equation  .  10 

Calibration  Procedure  . 11 

Data  Collection  Procedure  . 12 

Procedure  Validation . 15 

3—  Characterization  of  Soils  . 19 

Physical  Properties . 19 

Soil  Chemistry . 21 

Laboratory  procedures . 21 

Phase  composition  and  shape  parameters  . 23 

Bulk  chemical  composition . 28 

Chemistry  and  electrical  conductivity . 30 

4—  Analyses  of  Representative  Data . 31 

Dispersive  Behavior  . 31 

Sand . 41 

Silt . . 42 

Clay . 43 

A  hypothetical  sensor  design . 43 

Moisture  Effects . 44 

5—  Summary  of  Observations  and  Conclusions . 54 


References  . 

Appendix  A:  Data  Collection  and  Data  Presentation  Software . A1 

Appendix  B:  Electrical  Properties  of  Soils  Plotted  Against  Frequency  .  .  .  B1 

Appendix  C:  Electrical  Properties  of  Soils  Plotted  Against  Volumetric 

Moisture  Content . Cl 

SF  298 

List  of  Figures 


Figure  1 .  Dielectric  loss  mechanisms  for  heterogeneous  moist 

materials  . 7 

Figure  2.  Dielectric  property  measurement  apparatus  block  diagram  ....  9 

Figure  3.  Brass  square  cross-section  coaxial  sample  holders . 10 

Figure  4.  Filling  a  sample  holder  with  soil/water  mixture  . 13 

Figure  5.  Sample  holder  placed  in  sink/source  with  cables  attached  ...  14 

Figure  6.  Complex  dielectric  constant  for  air,  100-mm  holder  . 15 

Figure  7.  Complex  dielectric  constant  for  pure  water,  20  °C,  100-mm 

holder . 16 

Figure  8.  Complex  dielectric  constant  for  pure  water,  20  °C,  15-mm 

holder .  17 

Figure  9.  Soil  gradation  curves . 20 

Figure  10.  Electrical  properties  of  a  sand  . 32 

Figure  1 1 .  Electrical  properties  of  a  silt . 35 

Figure  12.  Electrical  properties  of  a  clay  . 38 

Figure  13.  Electrical  properties  of  all  soils  at  78  MHz . 45 

Figure  14.  A  bilinear  fit  to  permittivity  data  collected  at  78  MHz, 

20  °C  . 48 

Figure  15.  Least  square  polynomial  fits  to  permittivity  data  collected 

at  78  MHz,  20  “C  . 49 

Figure  16.  Permittivity  mixing  model  applied  to  a  sand,  78  MHz, 

20  °C  . 52 

Figure  17.  Permittivity  mixing  model  applied  to  a  silt,  78  MHz, 

20  °C  . 53 

Figure  18.  Permittivity  mixing  model  applied  to  a  clay,  78  MHz, 

20  °C  . 53 


List  of  Tables 


Table  1.  Physical  Characteristics  of  Soil  Samples .  19 

Table  2.  Measured  Concentrations  of  Mineral  Phases  With  Asso¬ 
ciated  Precision  Using  Quantitative  X-ray  Diffraction  Data  .  .  23 

Table  3.  Results  of  Chemical  Analyses  for  Individual  Particles  for 
Selected  Samples  Using  Standardless  EDX  Analysis  on 
SEM . 24 

Table  4.  Mineral  Phases  Assemblages  Given  in  Table  2  Adjusted  to 

100  Percent . 26 

Table  5.  Average  Sphericity  and  Roundness  for  Samples  Examined  in 

This  Study  . 27 

Table  6.  Results  of  Chemical  Analyses  for  Bulk  Samples  A 

Through  L  . 28 

Table  7.  Stoichiometries  Used  to  Calculate  Chemical  Compositions 

Based  on  Adjusted  Phase  Assemblages  From  Table  4 . 29 

Table  8.  Calculated  Chemical  Compositions  of  Samples  A  Through 
L  Based  on  Normalized  Quantitative  X-ray  Diffraction 
Data  and  Stoichiometries  for  Mineral  Phases  as  Indicated 
in  Table  7 . 29 


Preface 


The  report  herein  represents  one  in  a  series  of  reports  that  document  the 
results  of  a  basic  research  project  entitled  “In-Situ  Obstacle  Discrimination” 
(project  number  61102/AT22),  which  has  been  funded  by  the  U.S.  Army 
Corps  of  Engineers  Research  and  Development  Directorate  in  Washington, 
DC.  The  primary  goal  of  this  study  is  to  gain  a  fundamental  understanding  of 
the  mechanisms  associated  with  the  propagation  of  electromagnetic  energy 
through  soils  and  how  those  phenomena  affect  the  detectability  of  objects 
buried  beneath  the  surface. 

Dr.  John  O.  Curtis,  Environmental  Engineering  Division  (EED),  Environ¬ 
mental  Laboratory  (EL),  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES),  had  overall  responsibility  for  the  conduct  of  this  effort  and  was  the 
primary  author  for  this  report.  He  also  developed  the  experimental  procedure 
used  in  this  study  and  conducted  the  dielectric  data  analyses.  Dr.  Charles 
Weiss,  Concrete  Technology  Division,  Structures  Laboratory,  WES,  was 
responsible  for  the  petrographic  analyses  of  the  12  soils  and  made  use  of 
numerous  laboratory  and  analytical  tools  to  determine  the  chemical  composi¬ 
tions  of  the  soils.  Mr.  Joel  Everett,  EED,  conducted  all  of  the  dielectric 
property  measurements  and  created  the  initial  database  of  raw  results.  Soil 
samples  were  obtained  from  the  Soil  Test  Facility,  Geotechnical  Laboratory, 
WES,  where  many  of  the  physical  characteristics  were  also  identified.  Soil 
specific  surface  measurements  were  conducted  by  the  Quantachrome  Corpora¬ 
tion  in  Syosset,  NY. 

This  study  was  under  the  direct  supervision  of  Dr.  Raymond  L.  Mont¬ 
gomery,  Chief,  EED,  and  Dr.  John  Harrison,  Director,  EL. 
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Introduction 


Motivation  for  Study 


The  need  to  characterize  conditions  and  to  locate  anomalies  beneath  the 
terrain  surface  is  driven  by  numerous  scenarios  that  originate  in  both  the  mili¬ 
tary  and  private  sectors.  Soldiers  in  the  battlefield  need  to  locate  and  neutral¬ 
ize  buried  land  mines.  They  often  need  to  know  the  location  and  extent  of 
tunnel  systems  used  by  the  enemy.  Military  training  sites  require  frequent 
searches  for  unexploded  ordnance  (UXO)  to  ensure  the  safety  of  the  troops. 
Military  bases  that  are  scheduled  for  return  to  private  use  must  be  cleared  of 
UXO.  Chemical  contaminant  burial  or  dump  sites  at  both  military  and  non- 
military  sites  throughout  the  world  must  be  located  and  cleaned  up.  The 
extent  of  subsurface  contamination  at  burial  of  dump  sites,  or  even  storage 
sites,  where  the  containers  have  leaked  must  be  known  before  that  site  can  be 
remediated.  Civilian  construction  activities  can  be  made  safer  if  onsite  per- 
sonnel  have  an  accurate  picture  of  subsurface  conditions  including  utility  lines, 
cavities,  and  the  depth  to  the  water  table. 

One  of  the  most  effective  ways  of  safely  probing  beneath  the  terrain  sur¬ 
face  to  detect  anomalies  is  through  the  use  of  active  electromagnetic  systems 
such  as  ground-penetrating  radars  that  typically  operate  in  100-MHz  to  1-GHz 
frequency  range  (although  numerous  applications  exist  outside  of  this  range). 
Of  vital  importance,  then,  is  the  need  to  understand  how  the  electromagnetic 
waves  will  be  attenuated  through  the  subsurface  material  and  how  fast  they 
will  travel.  The  former  factor  determines,  when  coupled  with  the  sensitivity 
of  the  receiver  and  the  electrical  properties  of  material,  the  depth  at  which  an 
anomaly  can  be  detected.  The  latter  factor  is  particularly  important  for  accu¬ 
rately  determining  the  depth  to  the  anomaly  and  for  synthetic  aperture  process¬ 
ing  of  data,  where  keeping  track  of  the  phase  of  the  waves  is  of  paramount 
importance. 

Being  able  to  predict  the  performance  of  active  electromagnetic  systems  as 
tools  for  probing  subsurface  soil  conditions  requires  knowledge  of  the  com¬ 
plex  propagation  constants  of  the  soils  (to  be  discussed  in  the  following 
section)  and  how  they  are  affected  by  moisture,  temperature,  dry  density, 
chemical  composition  of  the  soils,  specific  surface,  grain-size  distributions, 
etc.  The  electrical  properties  of  the  soil- water-air-organics  mixture  define 
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both  the  phase  velocity  and  attenuation  information  that  is  required.  If  one 
can  develop  an  adequate  database  of  electrical  properties  for  a  broad  range  of 
soil  types,  moisture  conditions,  temperatures,  etc.,  researchers  and  designers 
can  make  informed  decisions  regarding  how  to  optimize  a  system  for  a  partic¬ 
ular  problem.  One  of  the  objectives  of  this  project  is  to  develop  such  a  data¬ 
base.  Hopefully,  these  data  will  complement  other  published  results.  While 
spotty  data  exist  at  high  frequencies,  there  have  been  thorough  studies  con¬ 
ducted  below  the  range  of  these  measurements  (Campbell  1990). 


Supporting  Studies 

The  volumetric  moisture  content  has  clearly  been  shown  to  be  one  of  the 
dominant  factors  in  controlling  the  electrical  properties  of  soils  over  a  broad 
range  of  frequencies  (Curtis  1993a,b,c).  Within  certain  frequency  ranges, 
other  important  elements  affecting  the  electrical  conductivity  of  soils  include 
the  mineralogic  composition,  the  type  of  clay  minerals  present,  and  possibly 
other  factors  such  as  sample  dry  density. 

An  excellent  set  of  soil  electromagnetic  property  measurements  were  con¬ 
ducted  at  the  United  States  Army  Electronics  Command  at  Fort  Monmouth, 
New  Jersey,  during  the  late  1960s  and  early  1970s  (Hulse,  Walker,  and 
Pearce  1969;  Walker,  Pearce,  and  Hulse  1969;  Walker  1970;  Hulse,  Walker, 
and  Pearce  1972;  and  Hulse  and  Walker  1972).  These  measurements  covered 
the  frequency  range  of  50  to  500  MHz  and  were  motivated  by  the  search  for 
ways  to  effectively  detect  enemy  tunnels  in  the  Vietnam  theater. 

One  of  the  conclusions  of  these  studies  was  that  soils  containing  high  con¬ 
centrations  of  smectite  group  minerals  were  observed  to  have  higher  conduc¬ 
tivities  than  those  not  containing  smectite  (or  expandable  clays),  primarily 
because  smectite  group  minerals  have  abundant  structural  water.  These 
researchers  even  attempted  to  develop  an  empirical  relationship  that  associates 
electrical  behavior  to  geochemical  composition.  For  example,  in  the  southern 
Atlantic  and  southeastern  areas  of  the  United  States,  the  conductivity  of  the 
soils  (ff)  at  frequencies  less  than  100  MHz  was  related  to  the  ratio  of  mont- 
morillonite,  a  smectite  group  mineral,  to  kaolinite  (m/k)  and  the  clastic,  frac¬ 
tion  by  weight  (c)  by  the  following  empirical  relationship  (Hulse,  Walker,  and 
Pearce  1972): 

ff  =  1  +  8.4  (m/A:)^-^  -  0.1c, 


Clastic  fraction  was  not  defined  by  the  authors  but  is  assumed  to  represent  the 
weight  fraction  of  the  soil  containing  mineral  and  rock  fragments.  The  units 
of  conductivity  are  mmhos/meter.  No  allowance  was  made  in  the  above  rela¬ 
tionship  for  moisture  content.  The  researchers  did  comment,  however,  that 
they  expected  each  clay  mineral  or  soil  system  to  have  its  own  unique  conduc¬ 
tivity  dependence  on  moisture  content.  Obviously,  determination  of  the 
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hydrous  phases,  especially  hydrated  clays,  is  going  to  be  important  to  the 
understanding  of  the  electrical  properties  of  soils. 

Another  excellent  source  of  both  electrical  and  magnetic  properties  of  soils 
is  the  United  States  Geological  Survey  in  Denver,  CO.  Since  the  late  1970s, 
numerous  measurements  made  there  and  by  colleagues  at  the  nearby  Colorado 
School  of  Mines  in  Golden,  CO,  have  been  conducted  as  part  of  an  overall 
effort  to  use  geophysical  methods  for  subsurface  characterization.  Some  of 
this  work  represents  summaries  of  measurement  programs  that  support  larger 
field  tests  (Olhoeft  and  Capron  1993).  Other  studies  reported  in  the  literature 
focus  on  the  physics  of  electromagnetic  propagation  through  soils  and  the 
chemistry  of  soil/water/organic/air  mixtures  that  produce  losses  in  the  propa¬ 
gating  signals  (Olhoeft  1985,  1987;  Olhoeft  and  King  1991;  Olhoeft  and 
Capron  1994).  The  results  of  these  studies  will  be  compared  with  the  results 
of  the  current  effort  in  a  later  chapter. 

Another  source  of  information  on  the  electromagnetic  properties  of  soils 
comes  from  research  conducted  at  the  University  of  California  at  Davis  since 
the  late  1960s.  Experiments  and  modeling  conducted  by  this  organization 
vary  from  electrical  analogues  used  in  an  empirical  fashion  to  duplicate  low- 
frequency  measurement  results  (Smith  1971)  to  sophisticated  analytical  models 
used  to  help  predict  geophysical  properties  of  soils  based  on  electromagnetic 
measurements  (Smith  and  Arulanandan  1981;  Arulanandan  1991,  1994). 

The  final  source  of  supporting  data  and  models  (although  certainly  not 
meaning  to  say  that  this  list  is  all  inclusive)  is  the  microwave  measurement 
program  led  by  Professor  Fawwaz  Ulaby,  first  at  the  University  of  Kansas 
during  the  early  1980s,  and  then  at  the  University  of  Michigan.  Although  a 
great  deal  of  the  work  conducted  by  this  group  focuses  on  radar  terrain  reflec¬ 
tances,  some  very  significant  dielectric  data  for  moist  soils  were  collected 
using  both  coaxial  devices  and  free-space  measurement  systems  (Hallikainen  et 
al.  1985).  An  attempt  was  also  made  by  this  group  to  develop  models  for 
predicting  dielectric  properties  that  were  both  theoretical  and  semiempirical 
(Dobson  et  al.  1985). 


Fundamental  Relationships 

The  relationships  that  follow  will  be  useful  in  interpreting  the  experimental 
results  as  will  be  discussed  in  a  later  chapter.  Supporting  references  include 
Jackson  (1975),  Stratton  (1941),  and  Bohren  and  Huffman  (1983). 

For  these  derivations,  plane  harmonic  wave  propagation  is  assumed  in  a 
lossy,  nonmagnetic,  unbounded  medium.  Therefore,  one  may  assume  that  the 
amplitude  function  may  be  written: 

giikx  -  ut) 
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where 


i  =  symbol  designating  an  imaginary  quantity  = 
k  =  P  +  ia  =  uN/c  (1) 

=  complex  propagation  constant 
(8  =  phase  constant 
a  =  amplitude  attenuation  factor 
03  =  radial  frequency 
N  =  complex  index  of  refraction 
c  =  velocity  of  light  in  a  vacuum 
X  =  a  space  coordinate 
t  =  time 
Furthermore, 

=  e  =  e'  -H  U"  (2) 


where,  as  stated  earlier,  e  is  the  relative  complex  dielectric  constant,  which, 
along  with  the  electrical  conductivity  from  Ohm’s  Law,  represents  the  electri¬ 
cal  properties  of  the  medium.  The  interpretation  of  these  properties  as  used  in 
this  study  is  that  the  conductivity,  a,  accounts  for  current  due  to  free  charged 
particle  motion,  while  the  imaginary  part  of  the  complex  dielectric  constant, 
e",  accounts  for  displacement  current  losses  (those  due  to  the  electric  polar¬ 
ization  of  the  medium).  When  both  conduction  and  displacement  currents  are 
considered,  one  finds  two  terms  in  Ampere’s  law  for  current  flow  that  repre¬ 
sent  losses  (or  a  shift  in  phase),  one  containing  the  electrical  conductivity  and 
one  containing  the  imaginary  part  of  the  dielectric  constant.  While  these  two 
terms  account  for  different  loss  mechanisms,  most  researchers  use  only  one 
term  or  the  other  to  identify  losses,  with  many  users  preferring  to  deal  with 
the  concept  of  electrical  conductivity.  In  MKS  units,  the  relationship  between 
the  two  quantities  is  taken  to  be 


where  the  units  of  conductivity  are  mhos/meter  (or  siemens/meter)  and  eo  is 
the  permittivity  of  free  space  (8.85  x  10'^^  farads/meter).  However,  in  the 
development  of  the  fundamental  relationships  used  to  interpret  laboratory  data. 
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it  was  found  more  convenient  to  account  for  the  dispersive  properties  of  the 
medium  completely  by  the  complex  dielectric  constant  (implicit  in  Equa¬ 
tion  2).  When  electrical  conductivity  is  desired,  Equation  3  will  be  utilized. 

Combining  the  square  of  Equation  1  with  Equation  2  and  equating  real  and 
imaginary  components  result  in  two  algebraic  equations  that  relate  the  ampli¬ 
tude  attenuation  factor  and  phase  constant  to  the  complex  dielectric  constant; 


c 


Solving  these  equations  for  the  amplitude  attenuation  factor  and  for  the  phase 
constant  results  in  the  following  expressions: 

1/2 

(4) 


(5) 


The  e"/e'  ratio  is  also  referred  to  as  the  loss  tangent.  Some  researchers  pre¬ 
fer  to  work  with  the  electrical  conductivity  (Equation  3)  in  place  of  the  dielec¬ 
tric  loss  term. 

Plane  waves  of  constant  phase  (but  not  necessarily  energy  or  signal)  will 
propagate  with  a  velocity 
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V  =  “  =  c 

/3 

e' 

2 

> 

r  1 2 
e" 

1 

(6) 


The  power  intensity  of  the  plane  electromagnetic  wave  decreases  exponen¬ 
tially  with  depth  of  penetration  by  the  factor,  e'^“*,  or,  in  one  unit  of  distance 
traveled,  a  decrease  of  e'^“.  Power  attenuation  expressed  in  decibels  per 
meter  can  then  be  written  as: 


PL  =  -8.6859  - 
c 

e' 

2 

> 

f  1  2 

f" 

1  +  ^  -  1 
[  e'  J 

1/2 


(7) 


Microwave  Loss  Mechanisms 


Before  examining  the  data  collected  in  this  study,  reviewing  several  of  the 
mechanisms  that  can  contribute  to  the  loss  in  energy  of  propagating  electro¬ 
magnetic  waves  would  be  helpful.  The  best  single  source  of  information  on 
this  topic  is  found  in  a  textbook  by  J.  B.  Hasted  (Hasted  1973).  In  particular, 
one  of  his  figures,  reproduced  here  in  Figure  1,  attempts  to  qualitatively 
identify  the  relative  magnitudes  and  approximate  spectral  range  for  several 
loss  mechanisms.  The  following  paragraphs  briefly  describe  those  mecha¬ 
nisms  that  apply  to  moist  soils  and  that  cover  the  range  of  frequencies  of  this 
study  and  may  prove  to  be  helpful  in  interpreting  much  of  the  data. 


Ionic  conductivity 

The  most  dominant  loss  mechanism  in  moist  heterogeneous  materials  at 
low  frequencies  is  that  due  to  ionic  conductivity,  the  a/(eoa))  term  in  Equa¬ 
tion  3.  Ionic  conductivity  results  when  salts  in  the  soil  are  dissolved  by  water 
and  a  path  is  available  through  the  pore  spaces  for  these  ions  to  travel  when 
subjected  to  an  external  electric  field.  Therefore,  low-frequency  ionic  conduc¬ 
tivity  in  moist  soils  has  to  be  related  in  some  way  to  the  degree  of  saturation 
and  how  the  pore  spaces  fill  with  water. 


Free  water  relaxation 

Individual  water  molecules  can  be  considered  to  be  electric  dipoles.  In  the 
presence  of  an  oscillating  electric  field,  the  dipoles  will  try  to  align  themselves 
with  the  field.  If  pure  water  was  placed  between  the  charged  surfaces  of  a 
parallel-plate  capacitor,  the  effect  would  be  an  enhancement  of  the  capacitance 
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C  ionic  conductivity 
DL  charged  doubie  iayer 
X  crystai  water  relaxation 
I  ice  relaxation 
MW  Maxwell-Wagner  effect 
S  surface  conductivity 
B  bound  water  relaxation 
W1  free  water  relaxation 


Figure  1 .  Dielectric  loss  mechanisms  for  heterogeneous  moist  materials 
(from  Hasted  1973) 


of  the  system  due  to  the  alignment  Of  these  dipoles  and  the  effective  reduction 
of  the  electric  field  (voltage  difference)  between  the  plates.  This  capacitive 
enhancement  is  measured  by  the  real  part  of  the  complex  dielectric  constant 
(also  referred  to  as  the  permittivity  of  the  material).  At  low  frequencies,  the 
permittivity  of  pure  water  is  about  80  times  that  of  air  or  a  vacuum). 

As  one  increases  the  firequency  of  the  oscillating  field,  the  dipoles  are 
forced  to  rotate  more  and  more  quickly.  At  some  frequency,  the  dipoles  will 
start  to  lag  behind  the  forcing  field  due  to  rotational  inertia  and  viscosity.  At 
this  point,  the  permittivity  of  the  substance  starts  to  decrease,  and  the  imagi¬ 
nary  part  of  the  dielectric  constant  starts  to  increase  to  reflect  the  increasing 
lag  between  the  electric  field  and  the  orientation  of  the  water  dipole.  As  the 
fi-equency  of  the  field  further  increases,  the  dipoles  get  to  the  point  where  they 
do  not  respond  to  the  field  at  all,  resulting  in  minimum  values  for  both  the 
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real  and  imaginary  parts  of  the  complex  dielectric  constant.  This  behavior  is 
referred  to  as  anomalous  dispersion  (“dispersion”  being  the  change  in  proper¬ 
ties  with  frequency).  It  is  also  referred  to  as  a  relaxation  phenomenon  with 
regard  to  the  measure  of  the  dipole’s  ability  to  return  to  its  equilibrium  posi¬ 
tion,  or  “relax,”  after  the  removal  of  an  electric  field. 


Bound  water  relaxation 

Experimental  evidence  exists  (Muir  1954;  Hoekstra  and  Delaney  1974)  that 
indicates  a  lowering  of  the  frequency  at  which  maximum  relaxation  losses 
occur  in  water  that  is  bound  to  solid  particle  surfaces.  Physically  this  makes 
sense  if  one  think.*;  of  the  water  to  solid  attraction  as  a  kind  of  added  viscosity . 
Hall  and  Rose  (1978)  proposed  bound  water  relaxation  as  the  dominant  loss 
mechanism  in  carefully  prepared  kaolinite  clays  in  the  10-  to  100-kHz  range 
of  the  spectrum;  however,  their  conclusions  seem  to  contradict  the  observa¬ 
tions  of  other  researchers. 


Maxwell-Wagner  effect 

Another  loss  mechanism  often  referred  to  in  the  literature  as  a  dominant 
mechanism  at  radio  frequencies  is  one  termed  the  Maxwell-Wagner  effect 
(Campbell  1988;  Bidadi,  Schroeder,  and  Pinnavaia  1988).  It  originates  in  the 
experimental  observation  that  there  is  dispersion  in  a  suspension  of  conducting 
particles  in  nonconducting  dielectric  media  much  like  that  observed  in  a 
parallel-plate  capacitor  filled  with  two  different  dielectrics.  The  Maxwell- 
Wagner  effect  is  seen  as  an  accumulation  of  charge  at  the  interface  of  dissimi¬ 
lar  materials  during  the  flow  of  current  through  heterogeneous  material  due  to 
a  discontinuity  in  dielectric  constant  values.  This  charge  buildup  is  time- 
dependent  and  results  in  dispersion  as  the  frequency  of  the  alternating  current 
increases  to  a  point  where  the  buildup  and  relaxation  cannot  keep  up. 


Study  Objectives 


The  objectives  of  this  study  are  to  describe  a  laboratory  procedure  for 
measuring  the  complex  dielectric  constant  of  a  material  (in  this  case,  moist 
soils)  and  to  report  the  results  of  those  measurements  on  a  large  number  of 
samples  in  various  forms  including  conductivity,  loss  tangent,  normalized 
phase  velocity,  and  power  loss  using  the  relationships  presented  above.  As 
will  be  seen  in  a  later  chapter,  these  parameters  are  a  strong  function  of  signal 
ft-equency,  moisture  content,  and,  in  some  instances,  soil  composition.  Sam¬ 
ple  temperature  (above  the  freezing  point  for  the  water  in  the  soil)  will  be 
found  to  be  a  much  weaker  contributor  to  the  dielectric  response  of  moist 
soils.  The  results  of  the  present  study  will  try  to  delineate  those  soils  in 
.which  the  greatest  attenuation  of  electromagnetic  power,  based  on  their  miner¬ 
alogy  and  chemical  composition,  will  occur. 
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2  Experimental  Technique  for 
Measuring  Complex 
Dielectric  Constants 


The  apparatus  used  to  collect  complex  dielectric  constant  data  is  a  modifi¬ 
cation  of  an  earlier  system  developed  at  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  (Curtis  1993b)  and  is  shown  in  block  diagram  form 
in  Figure  2.  The  system  is  built  around  a  Hewlett  Packard  8510C  Vector 
Network  Analyzer,  which  controls  a  synthesized  source  and  an  S-parameter 
test  set.  Specialized  flexible  cables  connect  the  test  set  to  unique  coaxial 
sample  holders  that  were  developed  at  WES.  The  holders  (shown  in  Figure  3) 
were  machined  from  brass  stock  and  have  a  cross  section  of  7.5  by  7.5  rnm. 
Sample  holders  were  manufactured  with  both  15-  and  100-mm-length  cavities 
to  ensure  that  valid  data  could  be  collected  at  the  extremes  of  the  frequency 
sweep.  Connection  to  the  test  set  cables  is  made  through  SMA  connectors  on 
the  ends  of  the  sample  holders  and  a  pair  of  3.5-mm  adapters.  The  complete 
apparatus  and  the  methodology  under  which  it  was  operated  are  somewhat 
similar  to  those  reported  in  another  recent  study  of  the  electrical  properties  of 
soils  (Wensink  1993). 


NETWORK  ANALYZER 

_ ^ 

CONTROLLER 

(HP8510C) 

^ - 

(HP9000) 

u  u 


Figure  2.  Dielectric  property  measurement  apparatus  block  diagram 
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Figure  3.  Brass  square  cross-section  coaxial  sample  holders 

Control  of  the  network  analyzer  measurement  system,  creation  of  raw  data 
files,  and  plotting  of  complex  dielectric  constants  are  accomplished  by  an 
HP-9000-362  instrument  controller.  A  copy  of  the  WES-developed  code, 
written  in  HP  BASIC,  can  be  found  in  Appendix  A.  Further  processing  and 
display  of  data  took  place  on  an  off-line  PC  using  the  SAS  software  system 
(SAS  Institute,  Inc.  1988). 


Governing  Equation 

The  synthesized  source  has  the  capability  of  generating  signals  that  span  a 
45-MHz  to  26.5-GHz  range.  Complex  reflection  and  transmission  coefficients 
(Sji  and  S21,  respectively)  are  measured  by  the  test  set.  Assuming  plane  wave 
behavior  within  nonmagnetic  samples,  one  can  show  that  the  relative  complex 
dielectric  constant  (e)  can  be  calculated  from  the  following  governing  equation 
(Curtis  1993b): 


COS  — ' 

c 


1  +  (52i)2  -  (5ii)2 


where  d  is  the  length  of  the  sample. 
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The  right-hand  side  of  this  equation  is  a  vector  in  complex  space  that 
rotates  about  the  origin  as  frequency  changes.  As  it  crosses  the  real  axis,  both 
on  the  positive  side  and  on  the  negative  side,  it  passes  through  or  near  branch 
cuts  for  the  inverse  cosine  function  (Abramowitz  and  Stegun  1972),  leading  to 
instabilities  in  the  calculation  of  the  argument  on  the  left.  As  will  be  seen 
later,  plots  of  the  real  part  of  the  complex  dielectric  constant  against  frequency 
take  on  the  characteristics  of  anomalous  dispersion  (Jackson  1975)  at  these 
points  where  the  complex  vector  approaches  the  real  axis. 

A  physical  description  of  what  is  happening  within  the  test  sample  at  these 
points  of  instability  may  best  be  demonstrated  when  one  attempts  to  measure 
the  response  of  the  empty  sample  holder;  i.e.,  a  measurement  of  the  complex 
dielectric  constant  of  air.  The  relative  constant  for  air  has  a  real  part  equal  to 
unity  and  an  imaginary  part  equal  to  zero.  One  can  then  rewrite  the  argument 
of  the  cosine  function  in  the  governing  equation  as  a  function  of  wavelength 
and  sample  cavity  length  by  setting  the  dielectric  constant  to  unity  and  substi¬ 
tuting  the  definition  of  frequency  in  terms  of  wave  velocity  and  wavelength 


At  the  same  time,  experimental  results  show  that  the  complex  vector  from  the 
right-hand  side  of  Equation  8  rotates  through  the  real  axis  on  both  sides  with  a 
value  of  nearly  unity.  Now  the  inverse  cosine  function  evaluated  at  unity  is 
2kTr,  and  the  function  evaluated  at  negative  unity  is  (tt  -I-  Ikir).  In  other 
words,  every  integer  multiple  of  tt  is  accounted  for.  Then  the  argument  of 
the  cosine  function  can  be  set  equal  to  an  integer  multiple,  n,  of  tt,  leading  to 
the  conclusion  that  as  the  complex  vector  crosses  the  real  axis, 

d  =  n{\/2) 

Therefore,  the  instability  in  the  calculation  of  complex  dielectric  properties 
takes  place  whenever  the  sample  contains  an  integer  multiple  of  half 
wavelengths. 


Calibration  Procedure 

Prior  to  collecting  data  on  moist  soil  samples,  a  full  2-port  calibration  was 
performed  using  an  HP85052B,  3.5-mm  Calibration  Kit.  Under  normal  con¬ 
ditions,  each  calibration  would  hold  for  several  hours  before  instrumentation 
drift  would  force  the  system  out  of  acceptable  limits.  The  limits  were  defined 
by  the  amplitude  and  phase  measurements  made  when  the  two  S-parameter 
test  cables  were  directly  connected.  A  calibration  was  considered  to  be  valid 
when  the  amplitude  did  not  vary  more  than  0. 1  db  over  the  entire  range  of 
frequencies,  and  when  the  phase  did  not  vary  more  than  1  deg  over  the  fre¬ 
quency  range. 
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As  described  elsewhere  (Curtis  1993b),  one  other  correction  had  to  be 
made  to  the  data,  even  after  a  successful  calibration,  before  correct  values  of 
the  dielectric  constant  could  be  computed.  Because  of  the  requirements  for 
adapters  between  the  sample  holder  SMA  connectors  and  the  S-parameter  test 
set  cables  and  because  of  the  existence  of  the  SMA  connectors,  themselves,  an 
electrical  delay  had  to  be  added  to  each  parameter  (Sn  Md  S21)  to  account  for 
the  additional  electrical  path  length  between  the  calibration  planes  at  the  ends 
of  the  cables  and  the  surface  of  the  sample.  For  these  measurements,  the 
necessary  electrical  delay  for  each  parameter  was  found  to  be 
270  picoseconds. 


Data  Collection  Procedure 


All  of  the  soils  used  in  this  study  were  organic-matter  free  and  had  been 
processed  for  ease  of  handling  (oven-dried,  pulverized,  and  screened  to 
remove  large  grains).  Sample  preparation  began  by  mixing  sufficient  deion¬ 
ized,  distilled  water  and  soil  in  a  plastic  bag  to  provide  a  mixture  that  was 
nearly  saturated.  The  bag  was  sealed  and  left  several  days  to  ensure  adequate 
time  for  adsorption  of  the  water  to  soil  particle  surfaces.  While  it  might  be 
useful  to  establish  a  target  gravimetric  moisture  content  (weight  of  water/ 
weight  of  dry  solid),  sample  preparation  procedure  would  not  ensure  that  the 
same  gravimetric  moisture  could  be  achieved  in  the  sample  holder.  Further¬ 
more,  volumetric  moisture  (volume  of  water/volume  of  sample),  which  is  the 
most  meaningful  measure  of  moisture  for  interpreting  the  dielectric  properties 
of  mixtures,  can  not  be  controlled  at  all  by  this  experimental  procedure. 

Once  the  soil/water  mixtures  had  been  properly  cured,  the  mixture  was 
“spooned”  into  both  the  15-mm  and  the  100-mm  sample  holders  (Figure  4). 
Removal  of  most,  if  not  all,  air  pockets  was  accomplished  by  tamping  the 
mixture  with  the  sharp  end  of  the  sample  preparation  tool  and  by  gently  strik¬ 
ing  the  sample  holder  against  a  flat  surface.  A  sharp  knife  edge  was  used  to 
flatten  the  sample  surface  and  clear  the  upper  edge  of  the  sample  holder  of 
particles.  The  covers  were  placed  on  the  sample  holders,  and  a  careful  mea¬ 
surement  of  the  combined  holder  and  sample  mass  was  made  (the  empty  sam¬ 
ple  holder  mass  was  known).  While  it  was  not  possible  to  ensure  that  both 
the  15-mm  and  the  100-mm  samples  possessed  the  same  volumetric  moisture, 
approximately  the  same  mechanical  effort  was  made  to  prepare  each  sample. 
Actual  differences  in  sample  dry  density  (weight  of  dry  solid/volume  of  sam¬ 
ple)  and  their  effect  on  dielectric  properties  will  be  discussed  in  a  later 
chapter. 

As  soon  as  the  soil/water  mixture  was  placed  in  the  two  holders  and 
weighed,  the  holders  were  placed  into  the  heat  sink/source  (Figure  5),  which 
had  been  brought  to  a  specified  temperature,  usually  10  °C.  The  cables  from 
the  S-parameter  test  set  were  then  attached,  and  the  cover  of  the  sink/source 
was  tightened.  A  Caron  Model  2065  Refrigerated  and  Heated  Circulating 
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Figure  4.  Filling  a  sample  holder  with  soil/water  mixture 


Bath  was  used  to  circulate  an  antifreeze  solution  through  the  sink/source. 
Once  any  required  temperature  was  achieved  in  the  Caron  bath,  the  sample 
was  allowed  to  reach  thermal  equilibrium  by  waiting  at  least  10  min  before  a 
measurement  was  made. 

Collection  of  complex  reflection  and  transmission  data  and  subsequent 
calculation  of  complex  dielectric  constants  at  each  desired  temperature  (10, 
20,  30,  and  40  °C)  were  achieved  by  interactively  executing  the  code  con¬ 
tained  in  Appendix  A.  The  entire  frequency  range  of  45  MHz  to  26.5  GHz 
was  divided  into  800  steps.  At  each  frequency  step,  the  and  S21  data  col¬ 
lected  represent  an  average  of  64  readings.  Following  a  frequency  sweep  for 
whichever  sample  holder  (either  the  15-mm  or  the  100-mm  holder)  was  con¬ 
nected  to  the  S-parameter  test  set,  the  cables  were  disconnected  and  moved  to 
the  other  holder  and  another  sweep  of  frequencies  was  executed.  After  mea¬ 
surements  were  made  at  one  temperature  setting,  a  new  sample  temperature 
was  selected  on  the  Caron  circulating  bath  control  panel,  and  another  set  of 
measurements  was  made  after  thermal  equilibrium  was  allowed  to  be 
achieved.  As  long  as  the  sample  had  reached  thermal  equilibrium,  the 
stepped-frequency  measurements  made  in  this  study  were  highly  repeatable. 

A  separate  data  file  was  created  for  each  sample  holder  and  at  each  temper¬ 
ature.  A  complete  set  of  measurements  on  a  pair  of  samples  at  all  four  tem¬ 
peratures  took  3  to  4  hr. 
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Figure  5.  Sample  holder  placed  in  sink/source  with  cables  attached 

If  only  one  moisture  content  condition  was  desired,  the  sample  holders 
were  removed  from  the  sink/source  as  soon  as  the  last  temperature  measure¬ 
ments  had  been  completed  and  were  weighed  again  to  determine  how  much,  if 
any,  moisture  has  been  lost  during  the  measurement  procedure.  Then  the 
samples  were  scraped  and  flushed  from  the  holders  into  separate  containers 
whose  empty  weight  was  known,  the  samples  were  dried  in  an  oven  for  24  hr, 
and  the  dry  solid  weight  of  each  sample  was  determined.  This  allows  for  an 
after-the-fact  determination  of  sample  volumetric  moisture,  gravimetric  mois¬ 
ture,  and  dry  sample  density,  the  relationship  among  these  parameters  being: 

volumetric  moisture  =  gravimetric  moisture  *  dry  soil  density 

For  soils  that  did  not  have  a  high  clay  mineral  content,  it  was  often  possi¬ 
ble  to  collect  data  at  different  moisture  contents  by  partially  drying  each  sam¬ 
ple  after  a  complete  set  of  temperature  measurements.  This  was  accomplished 
by  either  exposing  the  sample  to  the  air  for  a  period  of  several  hours  or  by 
heating  the  sample  holder  with  a  hair  dryer  to  evaporate  some  of  the  moisture. 
Assuming  that  one  has  the  same  volumetric  moisture  content  in  each  pair  of 
sample  holders  during  the  first  set  of  measurements,  and  knowing  the  volume 
of  each  sample  holder,  one  can  attempt  to  achieve  similar  moisture  contents  in 
the  two  holders  by  monitoring  the  weight  of  each  sample  holder  during  the 
drying  process.  This  technique  for  controlling  moisture  contents  without 
disturbing  the  samples  was  not  very  successful.  With  hindsight,  it  is  probably 
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better  to  make  up  different  moisture  content  samples  than  to  try  to  control  the 
drying  process.  Furthermore,  there  is  no  assurance  that  the  samples  possessed 
a  uniform  distribution  of  moisture  after  being  partially  dried,  and  the  impact 
of  nonuniform  moisture  distribution  on  measurement  results  is  very  difficult  to 
predict  (a  model  might  consist  of  both  series  and  parallel  capacitances  whose 
weighting  factors  are  a  function  of  how  the  sample  dries). 


Procedure  Validation 


As  an  example  of  the  results  generated  by  this  experimental  procedure, 
consider  the  measurement  results  shown  in  Figure  6.  These  data  are  for  an 
empty  100-mm  sample  holder,  thus  representing  a  measurement  of  the  com¬ 
plex  dielectric  constant  for  air.  Note  that  the  experimental  results  are  very 
close  to  the  theoretical  results  of  a  real  part  equal  to  unity  and  an  imaginary 
part  equal  to  zero.  Note  also  the  frequencies  at  which  the  calculation 
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Figure  6.  Complex  dielectric  constant  for  air,  100-mm  holder 
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instability  shows  up  as  anomalous  dispersions.  Calculations  verify  that  these 
unstable  points  occur  at  frequencies  where  an  integer  number  of  half  wave¬ 
lengths  exist  in  the  sample  holder. 

Similar  results  are  shown  in  Figures  7  and  8  for  measurements  on  distilled, 
deionized  water  at  20  °C,  Figure  7  being  for  a  sample  length  of  100  mm  and 
Figure  8  being  for  a  sample  length  of  15  mm.  Low-frequency  results  for  the 
long  sample  holder  match  published  data  on  fresh  water  (Ray  1972)  quite  well 
up  to  a  frequency  value  of  around  3  GHz.  Short  sample  holder  results  repre¬ 
sent  a  pretty  good  match  to  the  expected  behavior  of  water  at  20  °C  at  the 
highest  frequencies. 

These  data  demonstrate  several  points  that  need  to  be  made.  The  first  is 
that  the  shorter  sample  holder  is  required  to  obtain  data  at  the  highest  frequen¬ 
cies  for  very  moist  samples  due  to  the  attenuation  of  the  signal  through  the 
longer  sample  length.  However,  the  limit  of  system  sensitivity  (about  -90  db) 


Figure  7.  Complex  dielectric  constant  for  pure  water,  20  ®C,  100-mm  holder 
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for  the  transmission  parameter,  S21,  is  reached  even  for  a  15-mm-length  sam¬ 
ple  of  water. 

A  second  conclusion  that  can  be  drawn  from  these  data  is  that,  given  com¬ 
parable  sample  dry  densities  and  volumetric  moistures,  the  short  sample 
holder  gives  results  that  deviate  from  the  long  holder  measurements  at  low 
frequencies.  This  can  be  due  to  at  least  two  possible  factors.  One  is  that  the 
system  may  be  incapable  of  accurately  measuring  phase  shifts  for  lossy  mate¬ 
rials  in  the  short  holders;  however,  this  can  not  be  because  measurements  of 
the  empty  holder  are  quite  accurate.  The  other  is  that  during  the  calculation 
of  dielectric  constants,  one  encounters  from  Equation  7  a  constant  factor  of 
the  speed  of  light  divided  by  the  product  of  the  sample  length  and  test  fre¬ 
quency.  This  constant  must  be  squared.  At  the  low  end  of  the  frequency 
sweep  with  the  shorter  sample  holder,  this  factor  squared  has  a  value  of 
6,177.2;  whereas,  for  the  longer  sample  holder,  the  squared  factor  has  a  value 
of  119.7,  over  50  times  smaller.  It  is  entirely  possible  that  the  calculations  of 
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the  complex  dielectric  constant  are  made  more  sensitive  at  low  frequencies  by 
the  large  constant  factor  associated  with  the  short  sample  holder. 

Regardless  of  why  the  15-mm  sample  holders  have  difficulty  reproducing 
known  results  at  low  frequencies,  the  fact  remains  that  they  do.  As  a  result  of 
this  observation  and  the  fact  that  moist  sample  attenuation  will  very  likely  be 
too  high  to  obtain  useful  data  at  high  frequencies  from  the  100-mm  holder,  all 
of  the  complex  dielectric  data  reported  in  a  later  chapter  were  screened  in  the 
following  way.  All  data  reported  at  frequencies  below  3  GHz  were  collected 
from  the  100-mm  holder;  most  data  above  3  GHz  came  from  the  15-mm 
holder  because  of  the  unpredictable  attenuation  of  the  signal  at  these  frequen¬ 
cies  in  the  long  samples;  the  only  exceptions  were  for  measurements  made  on 
so-called  “dry”  soils  using  the  100-mm  holder.  While  an  effort  was  made  to 
produce  moist  soil  samples  of  comparable  volumetric  moisture  in  both  holder 
sizes,  the  truth  is  that  there  will  be  different  volumetric  moistures  associated 
with  high-  and  low-frequency  data.  It  is  also  possible  that  dry  soil  density  for 
each  sample  will  be  different,  and  it  may  be  possible  to  examine  the  impact  of 
that  variable. 

The  text  on  Figures  7  and  8  indicates  sample  lengths  of  97  mm  and 
14  mm,  respectively.  These  lengths  represent  a  final  calibration  adjustment 
through  the  calculation  of  the  complex  dielectric  constant  that  results  in  the 
best  match  to  accepted  values  for  water.  There  is  some  physical  basis  for 
using  sample  lengths  in  the  calculations  that  are  shorter  than  the  cavity  length 
in  that  the  SMA  connectors  have  dielectric  sleeves  around  the  center  conduc¬ 
tor  that  extend  a  short  distance  into  the  sample  cavity.  An  analytical  effort  to 
model  the  true  sample  holder  geometry  and  to  predict  an  effective  sample 
length  is  beyond  the  scope  of  this  project. 
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3  Characterization  of  Soils 


Physical  Properties 

Table  1  is  a  summary  of  the  physical  properties  of  12  different  soils  chosen 
for  this  study.  These  soils,  provided  by  the  Soil  Test  Facility  of  the  Geotech¬ 
nical  Laboratory  at  WES,  were  selected  to  span  the  gamut  of  grain-size  distri¬ 
butions  that  might  occur  in  nature  and  include  such  extremes  as  poorly  graded 
clean  sand  to  “fat”  clays.  Among  the  entries  found  in  Table  1  are  three  col¬ 
umns  referred  to  as  “Percent  Sand,”  “Percent  Silt,”  and  “Percent  Clay.” 

These  numbers  were  arrived  at  by  inspection  of  Figure  9,  which  contains  all 
of  the  gradation  curves  for  these  12  soils.  Sand  was  taken  to  be  the  percent 
(by  weight)  of  whatever  material  did  not  pass  the  number  200  sieve.  Clay 
content  represents  the  percent  of  material  whose  particles  exhibited  an  equiva¬ 
lent  size  of  less  than  2  /xm  (Mitchell  1993).  The  last  column  of  Table  1  is  a 
record  of  sample  specific  surfaces,  measured  by  Quantachrome  Corporation  in 
Syosset,  NY,  using  nitrogen  gas  sorption  techniques  (krypton  gas  was  used  for 
soil  “I”  because  of  its  extremely  low  specific  surface). 


Table  1 

Physical  Characteristics  of  Soil  Samples 

Sample 

Description  (Unified  Soil 
Classification  System) 

(by  weight) 

Liquid 

Limit 

Plastic 

Limit 

Specific 

Gravity 

Specific 

Surface 

m^/g 

Percent 

Sand 

Percent 

Silt 

Percent 

Clay 

A 

Clay  (CH),  Light  Gray 

2 

22 

76 

289 

39 

2.34 

25.97 

B 

Sand  (SP),  Light  Gray 

98 

2 

0 

2.62 

0.66 

C 

Silty  Sand  (SM),  Brown 

88 

8 

4 

-- 

-- 

2.64 

3.64 

D 

Silty  Sand,  (SM),  Reddish 

Brown 

77 

9 

14 

2.66 

15.01 

E 

Clayey  Silt  (ML),  Brown 

TR 

93 

7 

32 

24 

2.59 

13.89 

F 

Clay  (CH),  Gray 

1 

48 

51 

68 

21 

2.73 

33.34 

G 

Clayey  Sand  (SC),  Dark  Brown 

55 

32 

13 

38 

17 

2.66 

33.18 

H 

Clay  (CH),  Gray 

2 

64 

34 

61 

20 

2.74 

40.49 

I 

Sand  (SP),  White 

100 

0 

0 

- 

- 

2.67 

0.03 

J 

Silt  (ML),  White 

0 

46 

54 

47 

30 

2.61 

9.29 

K 

Clayey  Silt  (ML),  Brown;  Trace 
of  Sand 

4 

89 

■1 

27 

22 

2.71 

12.5 

L 

Sand  (SP),  Brown 

99 

1 

0 

- 

- 

2.66 

0.31 
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Figure  9.  Soil  gradation  curves 


Soil  Chemistry 


As  stated  in  the  Introduction,  soil  chemistry  has  to  be  one  of  the  primary 
factors  in  determining  the  electrical  behavior  of  the  soil/water/air  mixtures 
found  in  nature.  The  following  paragraphs  outline  a  number  of  soil  chemistry 
investigations  conducted  within  the  Petrography  and  Chemistry  Group  of  the 
Structures  Laboratory  at  WES.  The  objective  of  these  studies  was  to  identify 
sample  mineralogy,  chemical  composition,  and  particle  shape  characteristics 
that  might  influence  the  electrical  response  of  moistened  soils. 


Laboratory  procedures 

Qualitative  X-ray  diffraction  (XRD)  was  used  to  determine  the  mineralogy 
of  each  sample  and  also  to  estimate  the  quantity  of  each  mineral  present.  The 
equipment  used  in  this  analysis  was  a  Philips  APD1700  Version  4.0  Auto¬ 
mated  Powder  Diffractometer  system  driven  by  a  Digital  Equipment  Corpora¬ 
tion  (DEC)  VAX  3100  computer.  In  preparation  for  X-ray  diffraction  of  the 
bulk  samples  as  random  powder  mounts  (a  sample  with  random  crystal  orien¬ 
tations),  a  portion  of  each  sample  was  ground  in  a  mortar  and  pestle  to  pass 
through  a  45-/im  (No.  325)  mesh  sieve.  To  determine  if  clays  were  present,  a 
slurry  of  the  powder  with  water  was  made,  suspended  on  a  substrate,  and 
dried  overnight;  X-ray  diffraction  patterns  were  collected  on  these  oriented 
samples.  To  determine  if  the  clays  were  expandable.  X-ray  diffraction  pat¬ 
terns  were  collected  after  exposure  to  an  ethylene  glycol  atmosphere  overnight 
at  room  temperature. 

Random  powder  mounts  were  prepared  for  quantitative  X-ray  diffraction 
(QXRD)  analyses.  Five  standards  were  made  using  the  phases  found  in  the 
samples  and  an  internal  standard  (Ce02).  The  procedure  used  to  quantify  bulk 
powders  using  internal  standards  is  given  in  many  sources  (Klug  and  Alexan¬ 
der  1974;  Bertin  1975).  The  standards  were  used  to  obtain  calibration  coeffi¬ 
cients  for  the  intensity  of  a  reference  peak  for  each  mineral  as  a  function  of 
the  proportion  of  that  mineral.  These  coefficients  were  used  to  calculate 
mineral  concentrations.  Each  sample  was  analyzed  twice,  and  the  results  were 
averaged. 

Scanning  electron  microscopy  coupled  with  energy-dispersive  X-ray  (EDX) 
chemical  analysis  also  was  used  in  this  study  to  identify  any  additional  phases 
present  in  a  form  or  amount  not  detectible  by  X-ray  diffraction  analysis.  This 
was  accomplished  by  examining  the  morphology  and  determining  the  chemi¬ 
cal  composition  for  individual  particles.  The  equipment  used  was  a  Hitachi 
S2500  scanning  electron  microscope  (SEM)  with  an  EDX  windowless  detector 
linked  to  a  SUN  SPARCstation  IPC  running  the  IMIX  analysis  software 
(Princeton  Gamma-Tech,  Princeton,  NJ). 

Visual  examination  of  each  sample  was  done  to  verify  the  QXRD  results. 

In  this  examination,  a  portion  of  each  soil  was  placed  on  a  glass  slide, 
immersed  in  oil  (refractive  index,  n  =  1.544),  covered  with  a  cover  slip,  and 
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examined  using  a  Nikon  petrographic  microscope.  Verification  of  the  QXRD 
results  was  done  on  Samples  B,  C,  D,  G,  and  K.  This  procedure  consisted  of 
scanning  the  slide  sample  under  the  petrographic  microscope  along  parallel 
tracks  and  visually  identifying  the  type  of  mineral  grains  within  the  field  of 
view  until  a  total  of  200  grains  had  been  identified  and  tabulated.  Samples  A, 
E,  F,  and  H  could  not  be  point-counted  because  the  particle  grains  were 
agglomerations  of  various  minerals.  Samples  I  and  J  was  virtually  mono- 
mineralic,  and  these  samples  also  were  not  counted. 

Computer-based  image  analysis  was  used  to  characterize  the  size  and  shape 
parameters  of  each  soil.  In  this  analysis,  a  portion  of  each  sample  is  photo¬ 
graphed,  digitized,  and  analyzed  for  sphericity  and  roundness.  Particles  that 
were  touching  the  edge  of  the  image  or  touching  another  particle  were  not 
analyzed.  Organic  debris  such  as  seen  in  sample  C  also  were  ignored.  The 
equipment  used  in  this  analysis  included  a  Kodak  Eikonix  scanner  connected 
to  a  DEC  VAXstation  II/GPX  computer.  The  data  were  collected  and  ana¬ 
lyzed  using  the  Icometrics  version  3.0  software  package  (Everest  Technologies 
1991). 

Differential  scanning  calorimetric  measurements  were  made  to  ascertain  the 
presence  of  hydrous  phases.  The  equipment  used  was  a  Perkin-Elmer  DSC7 
differential  scanning  calorimeter  (DSC).  The  DSC  measures  the  amount  of 
energy  absorbed  or  evolved  by  the  sample.  In  these  analyses,  10  mg  of  each 
sample  was  placed  in  an  aluminum  pan,  covered  with  an  aluminum  lid,  and 
heated  from  20  to  200  °C  at  a  rate  of  10  °C/min. 

The  presence  of  the  acid-soluble  carbonate  minerals  calcite  and  dolomite 
was  determined  by  dissolution  in  HCl.  In  this  procedure,  5.0  g  of  sample  was 
put  into  100  mL  of  a  10-percent  solution  of  HCl  (by  volume)  and  then  stirred 
continuously  at  room  temperature  for  at  least  1  hr.  The  sample  was  then 
filtered,  dried,  and  weighed.  The  result  reported  for  this  procedure  is  the 
percent  weight  loss  subsequent  to  treatment.  The  resultant  sample  was  exam¬ 
ined  using  XRD  to  verify  that  the  carbonate  phases  had  been  removed. 

Bulk  chemical  composition  of  each  sample  was  determined  using  an 
energy-dispersive  X-ray  analyzer.  The  samples  were  prepared  by  making  a 
pellet  using  a  portion  of  the  <45-/im  fraction  and  a  methylcellulose  binder. 
The  equipment  used  was  a  Fisons-Kevex  771  EDX  chemical  analyzer  .for  bulk 
samples  running  the  XRF  Toolbox  II  version  4.4  software  on  a  Visions  486 
computer.  These  analyses  were  made  using  the  EXACT  analysis  procedure 
(Fisons-Kevex  1990)  and  a  plastic  clay  standard  (Source  Reference  Material 
from  NIST,  SRM  98b).  In  this  procedure,  the  standard  is  analyzed  and  the 
intensity  obtained  for  each  element  is  calibrated  to  the  known  chemical  com¬ 
position.  Subsequent  analysis  of  each  sample  is  given  relative  to  the  intensi¬ 
ties  for  the  standard. 
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Phase  composition  and  shape  parameters 


The  results  of  qualitative  X-ray  diffraction  were  used  to  select  standards 
for  quantitative  mineralogic  determinations.  All  phases  were  incorporated  in 
the  QXRD  analysis  with  the  exception  of  cristobalite,  which  was  present  in 
sample  A.  Due  to  the  lack  of  a  suitable  standard  of  cristobalite  the' quantity  of 
this  mineral  was  estimated  at  15  percent  of  the  total  based  on  the  peak  height 
of  the  bulk  pattern.  Results  of  QXRD  are  given  in  Table  2.  Several  factors 
could  contribute  to  the  totals  deviating  from  100  percent.  These  factors  could 
include  the  presence  of  additional  phases  or  poorly  crystalline  phases  not 
identified  by  XRD,  differences  in  extent  or  perfection  of  crystallinity  between 
minerals  in  the  samples  and  in  the  standards,  differences  in  total  chemical 
composition  between  samples  and  standards,  and  associated  error  in  the 
technique. 


Table  2 

Measured  Concentrations  of  Mineral  Phases  With  Associated  Precision  Using 
Quantitative  X-ray  Diffraction  Data 


Sample 


Quartz 


Smectite 


Calcite 


K-Feldspar 


Na-Plagioclase 


Mica 


Kaolinite 


Dolomite 


Total 


Cone.  Free.  Cone.  Free.  Cone.  Free.  Cone. 

%  abs  %  abs  %  abs  % 


Free.  Cone.  Free.  Cone.  Free, 

abs  %  abs  %  abs 


3.1  1.9 


25,6  2.2 


100.7  19.3  92.3  17.6  56.7  11.2  49.9  9.8  27.7 


Sample 


Quartz 


Smectite 


Calcite 


Cone.  Free.  Cone.  Free.  Cone.  Free.  Cone. 

%  abs  %  abs  %  abs  % 


Free.  Cone.  Free.  Cone.  Free, 

abs  %  abs  %  abs 


28.0  5.8  32.9  6.8  1 1 1 3.9  |  21 .3  |  1.0 


3.1  1.0  22.1  1.8 


0.9  40.8 


8.1  135.1  25.3 


K’Feldspar 

1.4 

0.4 

Na-Piagioclase 

3.0 

0.2 

2.6 

0.1 

Mica 

0.1 

0.0 

Kaolinite 

0.2 

1.3 

0.3 

Dolomite 


Total 


0.1  0.0 


114.0 


^  Cristobalite  also  was  present 
standard  could  not  be  obtained. 


in  this  sample,  but  was  omitted  from  the  quantitative  analysis  because  a  suitable 


Chapter  3  Characterization  of  Soils 


The  large  error  associated  with  the  quartz  measurement  could  explain  some 
of  these  discrepancies,  but  even  assuming  that  all  of  the  error  is  associated 
with  this  measurement,  the  totals  either  are  too  high  or  too  low.  Also,  it  is 
unlikely  that  an  error  in  the  calibration  constant  calculated  from  the  standards 
would  manifest  itself  in  overestimation  of  quartz  in  some  samples,  while 
underestimating  it  in  the  others. 

Because  the  method  used  for  QXRD  is  based  on  the  individual  peaks  of 
XRD  patterns  of  minerals,  variation  in  crystallinity  and  resulting  peak  heights 
will  be  manifested  in  the  quantitative  results.  For  example,  if  a  given  mineral 
sample  is  more  perfectly  crystalline  than  the  standard,  then  an  equal  amount 
of  that  mineral  in  the  sample  will  be  overestimated  in  the  analysis.  This  may 
explain  the  overestimation  of  quartz  in  some  of  the  samples. 

Differences  in  total  chemical  composition  between  the  minerals  in  the  sam¬ 
ples  and  those  of  the  standards  may  account  for  much  of  the  error  in  the 
analyses.  For  example,  the  smectite  standard  used  was  an  Fe-poor  Ca  mont- 
morillonite.  If  the  smectite  in  the  samples  either  contained  more  Fe  or  had 
some  substitution  of  Na  for  Ca  in  the  clay  interlayer,  this  would  lead  to  errors 
in  the  quantitative  analysis.  The  difficulties  associated  with  quantification  of 
clay  minerals  are  discussed  by  Zussman  (1967). 

To  verify  the  presence  of  Fe  and  a  variable  interlayer  cation  content  in  the 
structure  of  the  montmorillonite,  EDX  standardless  analysis  of  individual 
particles  was  undertaken.  The  results  of  these  analyses  are  given  in  Table  3. 
The  chemical  analyses  for  samples  A,  F,  G,  and  H  are  typical  for  those  sam¬ 
ples  and  have  similar  chemical  compositions  to  dioctahedral  smectites.  The 
analysis  for  sample  J  is  representative  of  a  kaolinite.  It  is  given  as  a  contrast 
to  the  smectite  analyses.  From  these  analyses,  it  is  obvious  that  the  smectites 
have  abundant  structural  Fe  and  a  variable  interlayer  cation  content.  These 
results  support  the  idea  that  the  smectites  were  underrepresented  in  the  QXRD 
analysis  because  the  clays  were  so  different  chemically. 


Table  3 

Results  of  Chemical  Analyses  for  Individual  Particles  for  Selected 
Samples  Using  Standardless  EDX  Analysis  on  SEM 

Composition 

A 

F 

G 

H 

J 

Si02 

68.3 

65.4 

59.0 

66.8 

55.2 

AI2O3 

23.1 

27.8 

30.8 

23.0 

43.8 

MgO 

3.8 

0.8 

1.1 

3.2 

0 

Na20 

2.8 

0 

0 

0 

0 

CaO 

0 

0.12 

0.7 

0.6 

0 

K2O 

0.2 

2.1 

1.4 

1.4 

0.1 

FejOg 

1.8 

3.5 

6.7 

4.7 

0 

Ti02 

0 

0.1 

0.4 

0.3 

0.9 
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Energy-dispersive  chemical  analyses  of  individual  smectite  particles  from 
samples  A,  E,  F,  G,  and  H  indicate  that  the  smectites  in  these  samples  have 
abundant  structural  Fe  and  variable  amounts  of  Na,  Ca,  and  K  in  their  struc¬ 
tures.  Compositional  differences  seem  to  explain  much  of  the  inconsistencies 
in  the  results,  as  those  samples  (A,  E,  F,  G,  and  H)  containing  smectite  have 
some  of  the  lowest  totals  in  Table  2.  Note  that  the  samples  that  did  not  con¬ 
tain  smectite  had  totals  either  close  to  or  greater  than  100  percent  (Table  2; 
Samples  B,  C,  I,  J,  and  L).  The  low  totals  in  Table  2  for  Samples  D  and  K, 
which  have  little  or  no  smectite,  appear  to  be  attributable  to  other  errors,  such 
as  differences  in  perfection  of  crystallinity.  No  additional  phases  were  indi¬ 
cated  by  observations  using  the  SEM. 

Results  of  the  quantitative  analyses  in  Table  2  were  adjusted  to  account  for 
the  discrepancies  discussed  previously.  These  adjusted  compositions  are  given 
in  Table  4.  For  those  samples  with  smectite  present  (A,  E,  F,  G,  and  H),  all 
of  the  difference  from  100  percent  was  attributed  to  smectite.  For  these  sam¬ 
ples,  no  other  phases  were  adjusted.  For  the  rest  of  the  samples,  the  phase 
compositions  were  normalized  to  100  percent.  For  samples  C  and  J,  which 
were  already  close  to  100  percent,  this  had  little  effect  on  the  amount  of  each 
mineral.  For  samples  D  and  K,  which  had  totals  further  from  100  percent, 
this  effectively  increased  the  amount  of  quartz  in  each  sample  by  26  to  32  per¬ 
cent,  but  had  smaller  effect  on  the  less  abundant  phases.  For  samples  B,  I, 
and  L,  which  had  totals  over  100  percent,  this  diminished  the  amount  of 
quartz  by  14  to  36  percent.  This  overestimation  of  quartz  in  the  original 
quantitative  analysis  indicates  that  the  quartz  standard  used  probably  had  a  less 
prefect  crystallinity  than  the  quartz  in  samples  B,  I,  and  L. 

The  visual  examination  of  these  samples  using  the  petrographic  microscope 
was  done  to  corroborate  the  XRD  results.  Examination  of  samples  I  and  L 
indicated  that  both  were  virtually  all  quartz.  In  agreement  with  the  XRD 
results,  a  trace  of  carbonate  was  observed  in  sample  I.  Sample  L  had  a  trace 
of  feldspar.  These  results  confirm  that  QXRD  overestimated  the  amount  of 
quartz  in  these  samples.  Furthermore,  the  results  of  point-counting  sample  B 
gave  52.5-percent  quartz,  45.5-percent  feldspar,  and  2-percent  mica.  This 
result  corroborates  the  amount  of  feldspar  as  determined  by  QXRD  and  that 
the  amount  of  quartz  is  overestimated. 

Visual  examination  of  samples  A,  E,  F,  G,  and  H  confirmed  the  presence 
of  clay  minerals.  The  clay  particles  in  these  samples  were  brownish  colored. 
The  type  of  clay  mineral  could  not  be  identified  by  this  examination.  As 
noted  previously,  samples  A,  E,  F,  and  H  could  not  be  counted,  but  one 
sample  containing  clay,  sample  G,  was  counted.  The  point-count  for  sample 
G  was  36-percent  quartz,  47-percent  clay,  3.5-percent  carbonate,  12.5-percent 
feldspar,  and  1-percent  mica  (percent  of  one  type  of  mineral  grain  out  of  total 
grains  counted).  These  results  compare  favorably  with  the  adjusted  QXRD 
results  (Table  4).  Although  all  the  samples  containing  clay  could  not  be 
counted,  this  one  result  confirms  that  the  QXRD  underestimated  the  amount  of 
expandable  clays.  Furthermore,  the  brownish  color  of  the  clays  indicates  the 
presence  of  Fe  in  the  structure. 


Chapter  3  Characterization  of  Soils 


25 


Table  4 

Mineral  Phases  Assemblages  Given  in  Table  2  Adjuste( 
100  Percent^ 

to 

Samite 

A 

B 

C 

D 

j _ 

F 

G 

H 

1 

J 

K 

L 

Quartz 

3 

71 

98 

89 

50 

28 

28 

33 

100 

1 

67 

99 

Smectite 

80 

0 

0 

0 

38 

67 

64 

62 

0 

0 

TR 

0 

Calcite 

TR2 

0 

0 

TR 

0 

0 

TR 

0 

0 

0 

3 

0 

K-Feldspar 

1 

23 

1 

6 

5 

TR 

5 

1 

0 

0 

3 

1 

Na-Plagioclase 

1 

6 

0 

0 

3 

0 

3 

2 

0 

0 

3 

0 

Mica 

0 

TR 

0 

0 

TR 

TR 

0 

TR 

0 

0 

TR 

0 

Kaolinite 

0 

0 

1 

5 

2 

5 

TR 

1 

0 

99 

2 

0 

Dolomite 

0 

0 

0 

0 

2 

0 

0 

0 

TR 

0 

21 

0 

Cristobalite 

15® 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

^  Sample  total  may  not  sum  to  1 00  percent  due  to  truncation  errors. 

^  TR  indicates  a  trace  amount  ( <  1  percent)  of  this  mineral  present. 

^  This  amount  is  estimated. 

Quantitative  X-ray  diffraction  analysis,  along  with  adjustments  described  in 
the  preceding  paragraphs,  provided  reasonable  phase  assemblages  for  samples 
in  this  study  (Table  4).  Using  these  data,  classification  of  the  soil  samples  can 
be  made  that  support  the  classifications  made  by  the  Soil  Test  Facility.  Sam¬ 
ple  A  is  a  smectite  bentonite  soil;  sample  B  is  a  quartzofeldspathic  soil;  sam¬ 
ples  C,  D,  and  L  are  quartzose  soils;  samples  E,  F,  G,  and  H  are  soils  rich  in 
quartz  and  smectite;  sample  I  is  virtually  a  pure  quartz  sand;  sample  J  is  a 
kaolinitic  soil,  and  sample  K  is  a  carbonaceous  quartzose  soil. 

Because  the  presence  of  hydrous  phases  is  an  important  factor  in  the  mea¬ 
sured  electrical  conductivity  of  soil,  differential  scanning  calorimetry  was  used 
to  verify  the  presence  of  these  phases.  Expandable  clays,  such  as  smectites, 
have  water  that  surround  (hydrate)  cations  within  their  structure,  which  is 
released  upon  heating  at  about  100  ®C  (Mackenzie  1957).  This  loss  of  water 
is  manifested  as  a  peak  due  to  the  endothermic  reaction.  The  DSC  result  for 
sample  A  clearly  shows  a  peak  at  97  °C,  which  is  ascribed  to  the  loss  of 
hydration  water.  The  peak  area  in  this  region  of  the  scan  is  directly  related  to 
the  amount  of  adsorbed  water  in  a  sample.  As  previously  noted,  samples  rich 
in  smectite  and  hydration  water  are  likely  to  have  higher  conductivities. 
Because  this  sample  had  the  largest  peak  due  to  loss  of  water  and  contained 
the  most  smectite,  it  probably  would  have  the  highest  conductivity  in  a  “dry” 
state.  Other  samples  (E,  F,  G,  and  H)  also  have  peaks  in  this  region  due  to 
loss  of  hydration  water.  These  samples  also  should  have  higher  electrical 
conductivities  than  the  rest  of  the  samples  devoid  of  hydration  water.  None  of 
the  other  samples  had  any  peaks  in  this  region.  These  results  agree  with  the 
XRD  results  that  smectite  is  present  in  samples  A,  E,  F,  G,  and  H. 
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Sample  J,  which  was  primarily  kaolinite,  did  not  disperse  properly  in  the 
oil  and,  therefore,  was  not  counted.  However,  as  already  discussed,  EDX 
analysis  on  the  SEM  indicated  that  the  chemical  composition  of  nearly  all  the 
particles  is  as  found  in  Table  3. 

For  samples  C,  D,  and  K,  there  was  good  agreement  between  the  point- 
count  results  and  the  adjusted  QXRD  results  (Table  4).  The  point-count  result 
for  sample  C  was  90-percent  quartz,  9-percent  feldspar,  0.5-percent  mica,  and 
1-percent  clay.  The  point-count  for  sample  D  was  87-percent  quartz, 

8-percent  clay,  3.5-percent  feldspar,  0.5-percent  mica,  and  1-percent  carbon¬ 
ate.  Again,  the  amount  of  carbonate  as  determined  by  dissolution  in  HCl 
(0.5  wt.  percent)  agrees  with  the  point-count  and  the  QXRD  results.  The 
point-count  for  sample  K  was  54-percent  quartz,  30.5-percent  carbonate, 
12-percent  feldspar,  1.5-percent  mica,  and  2-percent  clay.  Probably  due  to 
errors  in  counting,  the  amount  of  carbonate  by  point-count  was  slightly  higher 
than  the  adjusted  QXRD  value  (24  percent).  The  amount  of  carbonate  yielded 
during  the  dissolution  in  HCl  was  22.7  wt.  percent,  in  agreement  with  the 
adjusted  QXRD  results. 

Sphericity  and  roundness  were  calculated  from  data  obtained  using  image 
analysis  of  each  soil  (Table  5).  Roundness  relates  to  the  angularity  of  parti¬ 
cles.  When  the  curvature  or  angles  of  the  comers  of  a  solid  approach  that  of 
the  radius  of  the  inscribed  circle,  it  has  a  roundness  approaching  1 .  Spheri¬ 
city  states  quantitatively  how  nearly  equal  the  three  dimensions  of  a  particle 
are  (Folk  1968).  Samples  J  and  K  were  omitted  because  they  were  too  fine 
grained.  Values  for  sphericity  and  roundness  for  these  soils  range  from  0.4  to 
0.8  (Table  5).  These  parameters  were  calculated  from  the  data  obtained  from 
the  image  analysis  software.  The  sphericity  and  roundness  values  are  based 
on  the  analysis  of  particles  with  defined  sphericity  and  roundness  values 
(Krumbein  and  Sloss  1963).  Although  the  spread  in  values  for  sphericity  and 
roundness  is  small,  small  changes  in  these  values  have  a  great  efi^ect  on  the 
shape  characteristics  of  the  particles.  Sample  I,  which  is  virtually  all  quartz, 
has  the  largest  sphericity  (0.8)  and  roundness  (0.6)  values.  Copies  of  the 
images  used  for  this  analysis  are  on  file  and  available  on  request. 


Table  5 

Average  Sphericity  and  Soundness  for  Samples  Examined  in  This 
Study 

Sample 

A 

B 

C 

D 

E 

F 

G 

H 

1 

J 

K 

L 

Average  Sphericity 

0.6 

0.7 

0.6 

0.7 

0.6 

0.6 

0.7 

0.6 

0.8 

- 

0,7 

Average  Roundness 

0.4 

0.5 

0.4 

0.4 

0.4 

0.4 

0.5 

0.5 

0.6 

- 

- 

0.5 
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Bulk  chemical  composition 


Bulk  chemical  analyses  for  these  samples  are  given  in  Table  6.  Totals  for 
most  of  the  samples  are  close  to  100  percent.  The  low  totals  for  samples  E 
and  K  are  attributed  to  the  presence  of  dolomite  and/or  calcite  as  a  compo¬ 
nent.  The  low  totals  are  a  consequence  of  the  inability  of  the  energy- 
dispersive  chemical  analyzer  to  detect  the  presence  of  carbonate,  which  is  an 
important  component  in  these  samples. 


Table  6 

Results  of  Chemical  Analyses  for  Bulk  Samples  A  Through  L 

Composition 

A 

B 

C 

D 

E 

F 

G 

H 

■ 

J 

K 

SiO, 

74.3 

83.6 

91.9 

85.9 

73.8 

68.6 

75.4 

73.9 

80.3 

57.4 

60.2 

99.8 

19.7 

11.7 

5.5 

13.1 

11.5 

24.7 

15.5 

18.3 

2.6 

42.2 

9.6 

3.5 

MgO 

1.9 

0.3 

0,3 

0.9 

1.4 

0.6 

1.5 

0.0 

0.4 

2.2 

0.1 

mmm 

3.2 

3.2 

IQQ 

0.0 

1.1 

1.3 

1.4 

1.3 

0.3 

0.0 

1.5 

0.2 

0.6 

0.6 

0.3 

0.6 

0.1 

0.7 

0.4 

0.1 

0.3 

CaO 

1.3 

0.7 

0.1 

0.1 

1.0 

0.5 

0.6 

0.8 

0.1 

0.1 

0.1 

K2O 

0.3 

1.8 

0.1 

0.5 

1.0 

1.2 

1.4 

1.2 

0.0 

0.1 

0.2 

0.4 

0.8 

2.2 

3.3 

5.0 

3.5 

■a 

0.1 

0.5 

ID 

0.2 

Ti02 

HQ 

0.2 

0.6 

0.3 

KB 

0.8 

0.5 

0.7 

0.1 

1.9 

0.6 

0.2 

MnO 

0.0 

0.0 

0.0 

0.1 

0.1 

0.0 

0.1 

0.0 

0.0 

0.1 

0.0 

Total 

105.0 

102.4 

99.9 

102.6 

93.6 

104.1 

99.0 

102.7 

84.0 

102.6 

82.7 

104.5 

Although  sample  I  also  has  dolomite  as  a  constituent,  there  is  probably  not 
enough  present  to  account  for  the  disparity  observed.  This  interpretation  is 
corroborated  by  the  low  concentrations  of  CaO  and  MgO  (Table  6)  and  the 
loss  of  merely  0.26  wt.  percent  carbonate  by  dissolution  in  HCl.  The  X-ray 
diffraction  results  indicate  that  quartz  is  virtually  the  only  crystalline  phase  in 
this  sample,  and  the  chemical  composition  indicates  a  predominance  of  Si02. 
There  is  currently  no  satisfactory  explanation  for  the  low  compositional  total 
for  sample  I. 

The  hypothesis  that  the  normalized  phase  assemblages  from  Table  4  are 
consistent  with  the  bulk  chemical  compositions  was  tested  to  check  the 
assumptions  made  during  QXRD.  To  do  this,  the  normalized  phase  assem¬ 
blages  were  used  to  back-calculate  bulk  chemical  compositions.  A  chemical 
composition  for  each  sample  was  calculated  from  the  amount  of  each  phase 
given  in  Table  4  and  stoichiometries  for  the  mineral  phases  indicated  in 
Table  7.  The  results  of  these  calculations  are  shown  in  Table  8.  In  this  exer¬ 
cise,  the  amount  of  each  phase  with  a  concentration  less  than  1  percent  was 
assumed  to  be  0.  Comparing  data  in  Tables  6  and  8,  the  compositions  calcu¬ 
lated  from  the  adjusted  phase  assemblages  agree  reasonably  well  with  the 
results  of  the  bulk  chemical  analyses. 
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Table  7 

Stoichiometries  Used  to  Calculate  Chemical  Compositions  Based 
on  Adjusted  Phase  Assemblages  From  Table  4 

Quartz 

Si02 

Smectite 

^^0.25^®0.17ft^^2.88^®^^0.68^90.47J^^’7.71^l0.29^O20^^^Ul 

Calcite 

CaCOg 

K-Feldspar 

KAISiaOg 

Na-Plagiociase 

hl3o.8^^0.1^^^*3^8 

Mica 

K2Al4{Si6Al2)02o(OH)4 

Kaolinite 

AUSi40,o(OH)8 

Dolomite 

CaMg(C03)2 

Cristobalite 

Si02 

Table  8 

Calculated  Chemical  Compositions  of  Samples  A  Through  L  Based 
on  Normalized  Quantitative  X-ray  Diffraction  Data  (Table  3)  and 
Stoichiometries  for  Mineral  Phases  as  Indicated  in  Table  7 

Composition 

A 

B 

c 

D 

E 

F 

G 

H 

1 

J 

K 

L 

Si02 

68 

90 

99 

95 

80 

72 

73 

74 

100 

47 

72 

100 

AI2O3 

18 

5 

1 

3 

10 

16 

15 

15 

0 

39 

2 

0 

MgO 

2 

0 

0 

0 

1 

2 

2 

2 

0 

0 

5 

0 

Na20 

1 

1 

0 

0 

1 

0 

1 

1 

0 

0 

0 

0 

P2O5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CaO 

2 

0  . 

0 

0 

1 

1 

1 

1 

0 

0 

8 

0 

K2O 

0 

4 

0 

1 

1 

0 

1 

0 

0 

0 

1 

0 

Fe203 

6 

0 

0 

0 

3 

5 

5 

4 

0 

0 

0 

0 

Ti02 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

MnO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

HjO 

4 

0 

0 

1 

2 

4 

3 

3 

0 

14 

0 

0 

CO2 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

11 

0 

Total 

101 

100 

100 

100 

100 

100 

101 

100 

100 

100 

99 

100 

Measured  and  calculated  compositions  of  those  samples  with  smectite  in 
them  (A,  E,  F,  G,  and  H)  agree  closely.  The  good  agreement  between  the 
measured  and  calculated  chemical  compositions  for  sample  A  also  probably 
means  that  the  estimated  amount  of  cristobalite  was  not  excessive.  For 
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samples  B,  C,  D,  and  K,  the  predicted  amount  of  SiOj  in  each  sample  exceeds 
the  amount  measured.  It  is  probable  that  the  amount  of  quartz  was  overesti¬ 
mated  by  the  normalization  technique  used.  Other  normalization  schemes, 
such  as  keeping  the  amount  of  quartz  fixed  and  varying  the  other  phases,  are 
unsuitable  because  this  would  elevate  the  concentrations  of  the  other  cations 
above  the  amounts  actually  measured.  The  technique  used  is  imperfect,  but 
gives  a  justifiable  fit  between  predicted  and  measured  compositions. 

These  soil  chemistry  analyses  reveal  a  serious  problem  with  interpreting 
the  results  of  previous  studies  of  the  electrical  behavior  of  soils.  At  best, 
these  earlier  efforts  characterized  the  physical  properties  of  the  soils  being 
studied,  perhaps  listing  percentages  of  sand,  silt,  and  clay,  as  was  done  in 
Table  1.  Using  these  physical  properties  would  result  in  materials  E  and  J 
being  classified  as  silts  and  giving  no  indication  of  the  existence  and  character 
of  any  clay  minerals  that  may  be  present.  Table  4,  on  the  other  hand,  clearly 
shows  that  Soil  E  contains  a  significant  amount  of  expandable  clay  minerals, 
which  will  strongly  affect  its  electrical  behavior.  In  the  same  way.  Soil  J  is 
shown  to  be  pure  kaolinite,  a  nonexpandable  clay  mineral,  a  clay  mineral, 
nevertheless,  whose  surface  chemistry  is  quite  different  from  that  of  quartz, 
and  whose  electrical  behavior  is  likely  to  be  different  from  quartz. 


Chemistry  and  electrical  conductivity 

Because  the  electrical  conductivity  of  soils  is  strongly  controlled  by  the 
amount  of  water  and  because  the  amount  of  water  that  can  be  adsorbed  by  soil 
particles  is  greatest  for  smectite  clays,  estimates  of  the  electrical  conductivities 
of  these  soils  can  be  made  using  the  soil  chemistry  information  presented 
above.  Sample  A,  which  has  the  most  adsorbed  water  and  the  greatest 
amount  of  smectite,  should  provide  some  of  the  highest  low-frequency  con¬ 
ductivities  that  will  be  measured  using  the  apparatus  described  in  the  preced¬ 
ing  chapter.  The  other  samples  (E,  F,  G,  and  H)  with  observable  water  of 
hydration  and  abundant  smectite  concentration  also  should  provide  higher  elec¬ 
trical  conductivities  than  the  remainder  of  the  samples  whose  smectite  content 
is  low  or  nonexistent. 
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4  Analyses  of  Representative 
Data 


Well  over  500  measurements  of  complex  dielectric  constant  versus  fre¬ 
quency  were  collected  for  the  12  soils  described  in  the  previous  chapter.  Data 
were  recorded  using  samples  at  four  different  temperatures  (10,  20,  30,  and 
40  °C)  and  several  moisture  contents.  All  of  the  data  appear  in  the  Appen¬ 
dixes  B  and  C,  first  plotted  against  frequency  and  then  plotted  against  mois¬ 
ture  content  at  fixed  frequencies.  This  chapter  contains  samples  of  those  data 
that  will  be  examined  closely  in  an  attempt  to  draw  general  conclusions  about 
the  electrical  behavior  of  soils  as  a  function  of  physical  and  chemical  proper¬ 
ties  and  to  begin  to  explore  the  possibilities  for  modeling  such  behavior.  Data 
presented  in  this  section  are  for  measurements  made  at  a  single  temperature  of 
20  °C.  An  examination  of  all  of  the  data  reveal  only  a  small  effect  due  to 
changes  in  sample  temperature  over  this  entire  temperature  range,  an  obser¬ 
vation  that  has  already  been  reported  elsewhere  (Curtis  1993c). 


Dispersive  Behavior 


Examples  of  how  electrical  properties  of  these  soils  change  with  frequency 
are  shown  in  Figures  10-12,  which  contain  experimental  results  for  Soil  B 
(a  sand).  Soil  K  (a  silt),  and  Soil  H  (an  expandable  clay).  For  each  soil,  the 
figures  contain  plots  of  the  real  and  imaginary  parts  of  the  relative  complex 
dielectric  constant,  the  electrical  conductivity  of  the  soil/water/air  mixture 
calculated  from  Equation  3,  the  loss  tangent  (the  ratio  of  the  imaginary  to  the 
real  part  of  the  complex  dielectric  constant),  the  phase  velocity  of  the  plane 
wave  calculated  from  Equation  6  and  normalized  to  the  speed  of  light  in  a 
vacuum,  and,  finally,  the  loss  in  power  of  the  propagating  plane  wave 
expressed  in  decibels/meter  and  calculated  from  Equation  7 . 

As  mentioned  earlier,  very  little  soil  dispersion  data  exist  from  other 
sources  that  cover  a  broad  range  of  frequencies  and  that  are  supported  with 
adequate  physical  property  descriptions  to  allow  for  proper  comparisons. 
Nevertheless,  there  appears  to  be  good  correlation  between  the  data  collected 
for  this  study  and  that  reported  elsewhere  wherever  supporting  information  is 
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Figure  10.  Electrical  properties  of  a  sand  (Sheet  1  of  3) 
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Figure  10.  (Sheet  2  of  3) 
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Figure  1 1 .  (Sheet  2  of  3) 
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SOIL-K 


TEMPERATURE=20 


SOIL-H 


TEMPERATURE«20 


adequate  (Hulse,  Walker,  and  Pearce  1969;  Walker,  Pearce,  and  Hulse  1969; 
Delaney  and  Arcone  1982;  Hallikainen  et  al.  1985). 


Sand 

Consider  first  the  data  and  calculated  parameters  shown  on  Figure  10  for 
the  representative  sand.  The  first  thing  to  note  is  the  very  strong  dependence 
of  permittivity  (the  real  part  of  the  complex  dielectric  constant)  on  volumetric 
moisture  content.  The  greater  the  amount  of  water,  the  greater  is  the  number 
of  electric  dipoles  (water  molecules)  available  for  orientation  in  the  oscillating 
electric  field;  hence,  the  greater  is  the  permittivity.  This  strong  association 
with  moisture  with  be  examined  more  closely  in  a  following  section  dealing 
with  mixing  model  predictions  for  the  permittivity. 

Examination  of  the  plots  of  the  imaginary  part  of  the  dielectric  constant 
reveals  more  than  one  loss  mechanism  in  the  soil/water/air  mixtures  (Hasted 
1973).  As  discussed  in  the  Introduction,  low-frequency  losses  may  be  attrib¬ 
uted  to  conductivity  due  to  the  generation  of  ions  resulting  from  the  combina¬ 
tion  of  water  and  ionic  compounds  (salts)  in  the  soil.  Water  used  for  these 
studies  is  distilled  and  deionized.  Data  shown  in  the  next  section  reveals  that 
the  water,  by  itself,  does  not  exhibit  low-frequency  losses.  Salts  from  the  soil 
must  be  going  into  solution  with  the  added  water. 

Low-frequency  losses  may  also  be  attributed  to  the  Maxwell-Wagner 
effect,  in  which  effective  polarization  occurs  because  of  charge  buildup  on 
particle  surfaces  in  the  mixture  having  dissimilar  dielectric  properties.  Which 
mechanism  is  dominant  is  unclear  for  the  simple  reason  that  dielectric  prop¬ 
erty  measurements  were  not  made  at  frequencies  low  enough  to  determine 
whether  or  not  losses  are  anomalous  in  nature  (characteristic  of  Maxwell- 
Wagner  effects).  Anomalous  dispersion,  as  used  in  this  document,  refers  to  a 
decrease  in  permittivity  with  increasing  frequency  accompanied  by  a  maxima 
in  the  imaginary  part  of  the  dielectric  constant. 

Some  authors  have  effectively  argued  that  the  Maxwell-Wagner  effect  is 
dominant  in  soils  or  other  heterogeneous  mixtures  at  these  low  frequencies 
(Smith  1971;  Bidadi,  Schroeder,  and  Pirmavaia  1988).  In  a  recent  study  that 
made  use  of  laboratory  procedures  similar  to  the  ones  used  here,  one 
researcher  published  data  on  the  dielectric  response  of  natural  soils  between  1 
and  50  MHz  that  clearly  shows  a  Maxwell-Wagner-like  anomalous  dispersion 
although  he  never  referred  to  it  as  such  (Wensink  1993).  However,  if  an 
anomalous  dispersion  such  as  the  Maxwell-Wagner  effect  was  playing  a  part 
in  the  low-frequency  response  of  this  sand,  one  would  expect  to  see  a  decrease 
in  permittivity  along  with  a  decrease  in  the  imaginary  dielectric  constant  as 
frequency  increased.  This  trend  is  not  observed  in  the  sand. 

High-frequency  losses  are  undoubtably  due  to  the  classical  Debye  relax¬ 
ation  of  either  free  water  or  bound  water  molecule  dipoles  (Debye  1929; 
Hasted  1973).  Because  this  sand  is  relatively  free  of  clay  minerals,  its 
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high-frequency  response  is  probably  dominated  by  free  water  in  the 
soil/water/air  mixture. 

Although  the  imaginary  part  of  the  complex  dielectric  constant  is  associated 
with  a  shift  in  phase  of  a  propagating  electromagnetic  and,  hence,  is  related  to 
losses  in  the  mixture,  the  true  measure  of  loss  has  to  be  that  part  of  the  com¬ 
plex  propagation  constant  referred  to  as  the  amplitude  attenuation  factor 
(Equation  4).  Whereas  the  imaginary  part  of  the  dielectric  constant  indicates 
an  optimum  range  of  frequencies  (200  to  1,000  MHz)  in  which  losses  are 
minimized,  an  examination  of  the  power  attenuation  calculations  clearly  indi¬ 
cates  that  losses  increase  steadily  with  frequency.  This  is  because  the  increase 
in  frequency  is  the  dominant  term  in  Equation  7. 

Phase-velocity  calculations  provide  a  first  order  approximation  to  signal 
arrival  times.  This  would  be  very  important  in  an  iterative  scheme  for  pro¬ 
cessing  ground-penetrating  sensor  data  using  synthetic  aperture  techniques. 
Clearly,  phase  velocity  is  a  strong  function  of  the  permittivity  of  the  mixture. 


Silt 

Figure  1 1  contains  data  and  calculated  parameters  for  a  typical  silt  at  a 
sample  temperature  of  20  °C.  Two  significant  differences  between  silt  and 
sand  results  are  apparent  when  comparisons  are  made  at  similar  moisture 
contents.  The  first  is  that  the  permittivity  of  the  silt/water/air  mixture  is  not  a 
constant  over  a  broad  range  of  frequencies  as  was  true  for  the  sand.  This 
indicates  that  as  frequency  drops,  the  mixture  polarizes  more  easily.  The 
observed  behavior  is  compatible  with  an  anomalous  dispersion  model,  leading 
once  again  to  the  concept  of  the  Maxwell-Wagner  effect.  An  increased 
amount  of  charged  particles  are  available  in  this  silt  due  to  the  interaction  of 
water  and  the  dolomite  contained  in  the  soil.  It  would  be  the  movement  of 
these  ions  within  the  quartz  matrix  that  leads  to  charge  accumulation  at  parti¬ 
cle  boundaries  and  its  subsequent  relaxation. 

At  similar  volumetric  moisture  contents,  the  imaginary  part  of  the  dielec¬ 
tric  constant  is  greater  in  the  silt  than  in  the  sand  by  a  factor  of  about  five  at 
low  frequencies.  Although  the  data  are  not  as  complete,  there  appears  to  be 
little  difference  between  this  silt  and  sand  at  the  high-frequency  end  of  the 
spectrum.  This  is  probably  as  it  should  be,  because  the  mechanism  for  high- 
frequency  loss  is  dipole  relaxation  of  water  molecules.  If  the  amount  of  water 
is  basically  the  same,  and  if  there  is  no  reason  to  believe  that  the  water  mole¬ 
cules  are  being  constrained  in  the  silt  any  more  than  they  are  in  the  sand  (each 
soil  has  about  the  same  quantity  of  quartz),  then  the  behavior  of  the  two 
soil/water/air  mixtures  in  the  high-frequency  region  should  be  similar. 

Phase  velocities  in  the  silt  are  more  dependent  on  frequency  than  for  the 
sand.  Power  attenuation  in  the  silt  is  greater  at  low  frequencies  than  it  is  in 
the  sand  by  a  factor  of  about  five,  while  attenuation  in  the  two  soil/water/air 
mixtures  at  high  frequencies  is  about  the  same. 
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Clay 


Figure  12  shows  results  for  a  soil  that  contains  a  large  amount  of  expand¬ 
able  clay  minerals.  Although  there  is  only  a  limited  amount  of  high-frequency 
data  for  this  soil,  there  is  enough  along  with  the  low-frequency  results  to  draw 
several  conclusions.  The  first  is  that  at  low  frequencies,  similar  moisture 
content  samples  of  the  carbonaceous  silt  and  this  particular  clay  have  pretty 
much  the  same  complex  dielectric  response.  The  imaginary  part  of  the  con¬ 
stant  for  the  clay  is  about  twice  the  magnitude  of  that  for  the  silt  at  45  MHz. 
The  permittivity  of  the  clay  falls  off  more  quickly  with  increasing  frequency 
than  it  does  for  the  silt,  indicating  a  stronger  anomalous  dispersion  in  the  clay 
than  in  the  silt.  It  is  highly  probable  that  this  dispersion  is  enhanced  in  the 
clay  due  to  its  higher  specific  surface.  If  the  dispersion  mechanism  is  relax¬ 
ation  of  charges  built  up  on  solid  particle  surfaces,  and  there  is  more  surface 
area  for  charge  buildup  in  the  clay,  then  it  seems  reasonable  that  the  disper¬ 
sion  will  be  more  pronounced  in  the  clay. 

At  high-moisture  contents,  it  appears  that  all  three  representative  soils  have 
similar  dielectric  response.  This  could  be  explained  by  assuming  that  most  of 
the  water  in  the  soils  is  “free”  and  that  the  dielectric  response  of  the  mixture 
is  dominated  by  the  dipole  relaxation  of  the  “free”  water  molecules. 

At  low-moisture  contents,  there  is  clearly  a  difference  in  the  dielectric 
response  of  the  soils  due  to  how  the  water  and  soil  are  combined.  At  volu¬ 
metric  moistures  less  than  10  percent,  the  sand/ water/air  mixture  shows  the 
classical  dipole  relaxation  of  water;  hence,  in  the  sand,  most  of  the  water  must 
be  free  to  react  to  the  electric  field.  In  the  silt,  there  is  some  indication  of 
free-water  relaxation,  although  it  is  not  as  pronounced  as  in  the  sand.  As  for 
the  representative  clay,  a  measurement  made  for  what  was  thought  to  be  dry 
soil  (air  dried)  and  which  turned  out  to  be  about  10-percent  volumetric  water 
shows  a  low-frequency  loss  having  a  character  different  from  that  of  the  other 
soils  and  virtually  no  dipole  relaxation  at  high  frequencies.  The  water  in  the 
so-called  dry  clay  must  be  bound  to  the  clay  particles,  probably  within  the 
clay  interlayer  spacing  as  water  of  hydration  (Mitchell  1993),  but  not  so  much 
that  it  cannot  respond  to  the  applied  electric  field.  The  best  interpretation  of 
low-moisture  content  data  is  that  the  clay  exhibits  bound-water  relaxation 
(Hasted  1973). 

Results  of  calculations  of  the  effective  conductivity,  phase  velocity,  and 
power  attenuation  are  simply  reflections  of  the  dielectric  property  measure¬ 
ments.  Signal  attenuation  (power)  in  this  clay  will  be  about  two  times  the 
attenuation  of  the  silt  at  the  lowest  measurement  frequencies. 


A  hypothetical  sensor  design 

The  power  attenuation  curves  can  be  useful  as  a  tool  for  conducting  a  first 
order  design  of  a  ground-penetrating  electromagnetic  sensor.  For  example, 
the  sensitivity  of  a  receiver  is  such  that  it  can  detect  signals  (power  levels)  that 


Chapter  4  Analyses  of  Representative  Data 


43 


are  50  db  less  than  the  transmitted  signal.  Furthermore,  one  wishes  to  detect 
a  signal  from  a  buried  object  whose  reflected  signal  is  10  db  less  than  that  of 
the  incident  wave.  Finally,  the  sensor  operates  at  a  frequency  of  500  MHz 
and  suffers  no  loss  in  signal  power  at  the  air/ground  interface. 

The  calculated  power  attenuation  curves  for  the  representative  sand  and  the 
associated  range  of  moisture  levels  indicate  that  the  soil  could  attenuate  the 
signal  by  as  much  as  10  db  per  meter  of  travel  distance  for  a  volumetric  mois¬ 
ture  content  of  about  35  percent.  One  could  then  estimate  that  the  sensor 
would  be  able  to  detect  the  buried  object  at  depths  up  to  2  m  (greater  depths 
for  dryer  sand). 

Applying  the  same  parameters  to  the  design  of  a  sensor  that  might  operate 
in  the  representative  silt,  one  discovers  a  maximum  attenuation  of  about  30  db 
per  meter.  The  buried  object  discussed  above  would  not  be  found  below  a 
depth  of  about  two-thirds  of  a  meter  at  high  moisture  contents'. 

With  respect  to  the  representative  clay,  one  could  interpolate  the  calculated 
attenuation  curves  to  say  that  power  attenuation  at  500  MHz  and  about 
35-percent  volumetric  moisture  occurs  at  a  level  of  about  40  db  per  meter  of 
travel  distance.  The  end  result  of  the  design  calculation  is  that  the  particular 
buried  object  could  be  “seen”  in  the  clay  at  a  depth  of  no  more  than  one-half 
meter. 

Of  course,  these  example  designs  are  for  an  operating  frequency  of 
500  MHz,  somewhat  midrange  within  the  normal  operating  bands  of  ground- 
penetrating  radars.  Lower  frequency  systems  would  be  able  to  penetrate  much 
further  into  soils.  Although  the  data  from  which  these  values  of  attenuation 
were  calculated  were  collected  from  mixtures  of  soil  and  distilled,  deionized 
water,  which  has  a  low-frequency  behavior  much  different  from  fresh  water  or 
brine,  there  is  a  clear  indication  that  in  the  normal  range  of  operations  of 
ground-penetrating  radars  (<  1  GHz),  a  signal  should  be  able  to  penetrate 
almost  any  soil  to  a  shallow  depth  of  less  than  1  m  if  signal  strength  and 
sensitivity  are  both  adequate.  For  dry  sandy  soils,  measurable  signal  pene¬ 
tration  could  be  several  meters  or  more. 


Moisture  Effects 

From  the  representative  data  shown  in  the  previous  section,  one  immedi¬ 
ately  recognizes  the  very  strong  dependence  of  dielectric  properties  on  volu¬ 
metric  moisture  content.  Plotting  electrical  properties  versus  moisture  content 
at  some  frequency  of  interest  and  a  known  sample  temperature  will  provide 
further  insight  into  the  dielectric  response  of  soil/water/air  mixtures.  Con¬ 
sider  Figure  13,  which  contains  data  and  calculation  results  for  all  of  the  soils 
plotted  at  78  MHz  and  at  20  °C.  Although  low-frequency  electrical  data  are 
the  least  predictable  because  of  great  variation  in  the  loss  term,  several  useful 
observations  can  still  be  made. 
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Figure  13.  (Sheet  2  of  3) 
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Note,  first  of  all,  that  capacitive  properties  of  the  soils  (the  real  part  of  the 
dielectric  constant)  and  their  normalized  phase  velocities  hint  at  a  material- 
independent  response  when  considered  as  a  function  of  volumetric  moisture 
content.  Other  data  in  Appendix  C  for  higher  frequencies  show  an  even 
stronger  relationship  between  permittivity  and  volumetric  moisture  content. 
This  is  particularly  useful  if  one  considers  using  microwave  sensor  technology 
to  measure  moisture  in  the  soil. 

The  real  part  of  the  dielectric  constant  at  this  frequency  varies  nonlinearly 
with  moisture  just  as  has  been  reported  elsewhere  (Wang  and  Schmugge 
1980).  One  could  also  easily  imagine  a  bilinear  fit  to  these  data,  which  could 
be  interpreted  through  the  use  of  electrical  analogues  as  a  change  from  a 
model  of  series  capacitors  to  one  of  parallel  capacitors  at  some  critical  mois¬ 
ture  content  (Wang  and  Schmugge  1980;  Campbell  1988;  Curtis  1993c). 

These  models  reflect  the  creation  of  long-range  electrical  paths  in  the  soil/ 
water/air  mixture  at  the  critical  moisture  content  and  above.  For  the  78-MHz 
data,  a  bilinear  fit  (by  inspection)  would  result  in  a  critical  moisture  content  of 
about  24  percent  as  shown  Figure  14. 


Figure  14.  A  bilinear  fit  to  permittivity  data  collected  at  78  MHz,  20  °C 


One  could  easily  use  standard  statistical  techniques  to  perform  least-squares 
fits  to  such  data  as  is  shown  in  Figure  15,  where  one  second-order  polynomial 
fit  of  all  data  at  20  °C  is  compared  with  a  similar  fit  of  data  having  less  than 
36-percent  moisture  by  volume.  Because  natural  soils  normally  do  not  exhibit 
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Figure  1 5.  Least  square  polynomial  fits  to  permittivity  data  collected  at 
78  MHz,  20  °C 


volumetric  moistures  higher  than  this  level,  the  latter  fit  would  be  a  very  good 
representation  of  the  permittivity-volumetric  moisture  relationship  for  soil/ 
water/air  mixtures  at  78  MHz  and  20  °C. 

Another  option  for  fitting  a  model  to  the  permittivity-moisture  content  data 
in  Figure  13  is  an  exponential  relationship 

/  =J  gk*MC 

^  mixture  dry  soil 

where  MC  is  the  moisture  content  expressed  as  a  decimal.  Within  these  con¬ 
straints,  the  constant  factor,  k,  is 

d 

^  Ytater 

77 

^  dry  soil 

which,  for  the  data  shown  in  Figure  13,  have  a  value  of  about  3.45.  Check¬ 
ing  the  midpoint  of  the  exponential  model  fit,  one  sees  that  the  model  predicts 
a  permittivity  of  about  14,  whereas  the  data  show  an  average  value  on  the 
order  of  35.  An  exponential  model  is  not  adequate  for  the  composite  data. 
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Only  permittivity  data  collected  at  very  low  frequencies  and  high-moisture 
contents  exhibited  any  sort  of  variance  about  the  means  at  those  moistures. 
Note  also  that  all  of  the  soils  whose  permittivities  did  not  exhibit  predictable 
behavior  were  soils  that  contained  a  reasonably  large  amount  of  the  expand¬ 
able  clay  mineral,  smectite.  Because  permittivity  is  a  result  of  polarization 
phenomena  within  the  soil/water/air  mixture,  something  is  different  about 
what  causes  polarization  in  the  smectitic  soils  as  opposed  to  other  soil  types. 

In  sandy  soils  with  little  or  no  clay  content,  polarization  is  controlled  by  rota¬ 
tion  of  the  “free”  water  dipoles.  As  was  argued  in  the  previous  section,  wet 
silts  and  clays  are  both  likely  to  contain  ions  that  effectively  produce  polariza¬ 
tion  by  charge  accumulation  on  soil  particle  surfaces.  It  was  further  argued 
that  some  other  mechanism  had  to  be  responsible  for  the  enhanced  permittivity 
in  the  smectitic  soils,  clearly  evident  in  the  measurements  of  so-called  dry 
soil.  All  sources  seem  to  point  to  the  response  of  water  molecules  that  are 
bound  to  the  clay  mineral  particles,  and  whose  resonant  frequency  (inverse  of 
relaxation  time)  is  on  the  order  of  100  MHz  (Hasted  1973). 

Another  useful  observation  from  the  data  shown  in  Figure  13  can  be  made 
from  the  power  attenuation  plot.  About  an  order  of  magnitude  difference  in 
power  attenuation  at  any  moisture  content  over  the  range  of  soil  types  as  well 
as  an  expected  increase  in  attenuation  with  increasing  moisture  is  first  of  all 
seen.  The  upper  bound  of  these  data  could  be  used  as  a  conservative  design 
parameter  for  a  given  ground-penetrating  radar  system,  from  which  one  might 
conclude  that  power  loss  in  nearly  dry  soils  is  on  the  order  of  a  few  decibels/ 
meter,  while  saturated  soils  may  result  in  40  to  50  db/meter  of  loss.  Data  in 
the  appendixes  indicate  even  higher  losses  at  some  frequencies. 

Because  of  the  strong  correlation  between  soil/water/air  mixture  permittiv¬ 
ities  and  volumetric  moisture  demonstrated  by  the  composite  data,  one  might 
be  strongly  tempted  to  utilize  a  permittivity-mixing  model  based  on  electrical 
analogues  for  the  physical  mixtures.  In  other  words,  one  could  equate  the 
permittivity  of  each  component  in  the  mixture  to  the  capacitance  of  a  simple 
parallel  plate  capacitor.  These  capacitors  could  be  arranged  in  a  series  circuit, 
a  parallel  circuit,  or  in  a  more  cothplex  configuration.  The  effective  capac¬ 
itance  of  the  combined  elements  would  represent  the  permittivity  of  the  mix¬ 
ture.  The  most  general  form  of  such  a  mixture  model,  as  it  applies  to 
mixtures  of  dry  soil,  pure  water,  and  air,  is  (Curtis,  In  Preparation): 


”  ^soi^  soil 


6  e*  + 

water  water 


6  €* 
air  air 


where  the  subscripted  deltas  represent  volume  fractions  of  the  components  and 
the  primed  notation  for  the  real  part  of  the  dielectric  constant  has  been 
dropped.  One  implicit  assumption  to  this  kind  of  model  is  that  each  of  the 
constituent  materials  do  not  lose  their  individual  capacitive  properties.  In 
other  words,  there  can  be  no  chemical  interaction  among  the  mixture  compo¬ 
nents  that  changes  their  respective  dielectric  responses.  The  mixing  law  expo¬ 
nent,  k,  takes  on  the  special  values  of  + 1  if  the  electrical  paths  through  the 
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mixture  components  are  parallel  to  each  other  and  -1  if  the  current  is  forced  to 
pass  through  each  component  in  a  series  arrangement. 

Such  a  model  cannot,  however,  be  applied  to  the  composite  data  shown  in 
Figure  13  for  the  following  reasons.  First  of  all,  if  one  wished  to  fit  a  single 
mixing  model  to  the  composite  data,  the  volume  fraction  for.  water  would  have 
to  vary  from  0.0  to  1.0.  What  meaning,  then,  could  be  given  to  the  volume 
fraction  of  the  dry  soil  element  or  to  that  for  air?  Most  of  the  samples  from 
which  the  composite  data  are  derived  had  different  solid  fractions.  In  fact, 
high  moisture  content  in  these  data  is  synonymous  with  low  solid  fractions. 
Even  if  some  kind  of  average  solid  fraction  was  chosen  for  all  soil  types,  the 
model  would  still  fail  at  high-moisture  contents  because  of  the  physical  con¬ 
straint  that  the  sum  of  the  volume  fractions  must  be  equal  to  unity. 

The  only  sensible  way  to  use  a  mixing  model  like  that  above  is  to  apply  it 
to  a  series  of  mixtures  whose  solid  fraction  remains  fixed  and  whose  water 
fraction  is  allowed  to  vary.  The  maximum  allowable  moisture  fraction  would 
then  be  limited  to  the  conditions  for  saturation. 

Fortunately,  a  number  of  the  data  points  collected  for  this  study  do  repre¬ 
sent  measurements  made  of  samples  that  were  first  tested  at  near-saturation 
conditions  and  then  dried  and  retested  without  disturbing  the  physical  sample. 
In  fact,  each  of  the  soils  whose  dispersion  curves  were  shown  in  the  previous 
section  were  examined  this  way.  It  should  be  possible,  then,  to  test  the 
hypothesis  that  a  dielectric  mixing  model  such  as  that  shown  above  is  useful 
for  predicting  the  permittivity  of  mixtures  at  other  moisture  contents. 

The  application  of  the  generalized  mixing  model  to  establish  the  mixing 
law  exponents  for  each  soil  type  requires  that  one  additional  factor  be  taken 
into  account.  For  so-called  dry  soil  samples,  the  permittivity  measurement  is 
really  a  measurement  of  a  mixture  containing  dry  solid  material  and  air.  The 
model  must  apply  for  this  condition  as  well.  What  is  required,  then,  is  that  a 
value  for  permittivity  for  the  dry  solid  component  be  chosen  to  force  the 
model  to  generate  the  correct  “dry”  soil  measurement  result.  Actually,  this 
solid  fraction  permittivity  is  calculated  for  each  value  of  the  mixing  exponent 
using  the  measured  dry  solid  fraction,  an  air  volume  fraction  that  is  one  minus 
the  dry  solid  fraction,  and  an  appropriate  dry  soil  mixture  permittivity.  The 
results  of  calculating  solid  fraction  permittivities  using  this  technique  are  very 
much  in  line  with  published  data  on  rocks  and  compressed  soils  (Nelson, 
Lindroth,  and  Blake  1989;  Ulaby  et  al.  1990). 

Figures  16,  17,  and  18  show  the  results  of  fitting  the  generalized  mixing 
model  at  78  MHz  and  20  °C  to  the  volume  fraction  parameters  associated 
with  the  measurements  of  soils  B,  K,  and  H,  respectively.  The  permittivity  of 
pure  water  was  taken  to  be  78.96;  solid  fractions  were  0.627  for  soil  B,  0.531 
for  soil  K,  and  0.32  for  soil  H.  Superimposed  on  the  calculated  mixture 
model  responses  for  various  values  of  the  mixing  exponent  are  the  measured 
data  for  the  soil/water/air  mixmres. 
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Figure  16.  Permittivity  mixing  model  applied  to  a  sand,  78  MHz,  20  °C 


Conclusions  that  can  be  drawn  from  Figures  16-18  include  the  following 
statements.  First  of  all,  the  mixing  model  exponents  are  clearly  material 
dependent.  One  could  argue  that  this  observation  is  misleading  because  sam¬ 
ples  were  measured  at  different  dry  densities.  Unfortunately,  there  was  no 
control  over  this  variable,  and  the  concern  cannot  be  addressed. 

Another  conclusion  is  that  the  only  decent  correlation  between  data  and 
models  was  for  the  sandy  soil,  B.  Neither  the  silt  data  nor  the  expandable 
clay  soil  data  were  well  represented  by  a  single  mixing  exponent.  One  possi¬ 
ble  explanation  is  the  possibility  of  chemical  reactions  within  the  silt  and  clay 
due  to  the  addition  of  water.  As  noted  in  an  earlier  section,  chemical  reac¬ 
tions  could  account  for  unexpected  low-frequency  permittivity  and  conductiv¬ 
ity  results,  and  would  be  in  direct  conflict  with  the  model  assumption  that 
each  component  of  the  mixture  retain  its  unique  dielectric  properties. 

Because  of  the  extreme  variability  of  dry  soil  densities  in  the  real  world 
along  with  the  fact  that  water  and  some  soil  minerals  combine  to  produce 
unpredictable  dielectric  responses,  it  seems  unlikely  that  any  advantage  can  be 
realized  from  trying  to  fit  a  generalized  mixing  model  to  permittivity-moisture 
content  data. 
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5  Summary  of  Observations 
and  Conclusions 


Broadband  (45  MHz  to  26.5  GHz)  complex  dielectric  property  measure¬ 
ments  were  made  for  12  different  soils  having  a  wide  spectrum  of  physical 
and  chemical  properties.  Data  were  collected  at  numerous  volumetric  mois¬ 
ture  contents  and  at  four  distinct  temperatures  (10,  20,  30,  and  40  °C).  Pre¬ 
sentation  of  the  data  was  made  in  various  forms  including  not  only  values  of 
the  measured  complex  dielectric  constant,  but  also  computed  values  of  con¬ 
ductivity,  loss  tangent,  normalized  phase  velocity,  and  power  attenuation  for 
idealized  plane  wave  propagation.  These  data  were  plotted  against  frequency 
for  each  temperature  and  each  sample  volumetric  moisture  and  against  mois¬ 
ture  content  at  fixed  frequencies,  and  correlated  well  with  limited  data  avail¬ 
able  elsewhere  in  the  literature. 

Dispersion  data  revealed  separate  low-frequency  and  high-frequency  loss 
mechanisms  (a  rough  division  between  low-  and  high-frequency  behavior 
being  about  1  GHz).  High-frequency  losses  are  attributed  to  free-water  mole¬ 
cule  dipole  relaxation  and  appear  to  be  independent  of  soil  type.  Low- 
frequency  losses  seem  to  be  associated  with  the  movement  of  ions  in  the 
soil/water/air  mixtures  that,  due  to  the  presence  of  a  solid  matrix,  result  in 
charge  accumulation  on  particle  surfaces  and  an  effective  polarization  and 
relaxation.  Measurements  were  not  made  at  frequencies  sufficiently  low  to 
separate  the  contribution  of  this  effect  (known  as  Maxwell-Wagner)  from  pure 
ionic  conduction.  The  soils  containing  expandable  clay  minerals  demonstrate 
another  loss  mechanism  at  low  frequencies;  namely,  bound  water  relaxation 
due  to  the  presence  of  water  of  hydration.  Low-frequency  losses  are  defi¬ 
nitely  related  to  the  chemistry  of  the  soil  through  either  the  generation  of  ions 
by  salts  within  the  soil  or  by  the  presence  of  water  within  the  soil  mineral 
crystals  due  to  hydration.  Physical  properties  (dry  density,  moisture  content, 
grain-size  distributions,  specific  surface,  etc.),  alone,  are  not  sufficient  to 
predict  the  electrical  response  of  all  soils. 

The  permittivity  of  soil/water/air  mixtures  is  independent  of  soil  type  at 
■higher  frequencies  when  considered  as  a  function  of  volumetric  moisture 
content.  Even  at  the  lower  frequencies  examined  in  this  study,  a  good  fit  to 
the  permittivity  data  can  be  made  by  a  bilinear  function  with  a  break  at  about 
24-percent  volumetric  moisture.  Very  effective  models  for  predicting 
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permittivity  can  also  be  generated  by  least-square  polynomial  fits  over  a  rea¬ 
sonable  range  of  moisture  values.  Exponential  models  do  not  predict  permit¬ 
tivity  if  fitted  to  the  entire  range  of  moistures. 

Electrical  analogues  can  be  used  to  generate  permittivity  mixture  models 
for  data  in  which  the  moisture  content  was  allowed  to  vary  for  a  fixed  volume 
fraction  of  solid  particles.  Arguments  are  given  regarding  why  these  models 
cannot  be  applied  to  the  composite  data  for  all  soil  types.  The  mixture  models 
appear  to  be  effective  only  for  soils  in  which  chemical  reactions  (ionization, 
effective  bonding  of  water  by  various  mechanisms)  have  not  taken  place. 
Electrical  analogues  must  be  used  very  cautiously  when  attempting  to  model 
the  dielectric  response  of  moist  soils. 

As  with  any  useful  study,  attempting  to  answer  one  set  of  questions  gener¬ 
ates  another.  More  experiments  are  needed  to  better  understand  the  effect  of 
variable  minerology,  chemical  composition,  particle  size,  and  water  content  on 
the  dielectric  properties  of  soils.  The  following  suggestions  are  made  for 
future  studies. 

a.  To  assess  the  effect  of  particle  size  as  well  as  to  address  questions  of 
varying  sample  dry  density,  a  number  of  samples  of  the  same  soil  type 
could  be  prepared  for  testing  by  selectively  screening  the  soil  to  pro¬ 
duce  samples  with  different,  but  nearly  uniform,  particle  sizes. 

b.  To  address  the  structural  and  chemical  effects  of  expandable  clays  on 
the  permittivity  of  the  soil/water/air  mixture,  soil  samples  could  be 
prepared  so  that  each  contains  a  different  smectitic  clay  mineral. 

c.  To  eliminate  the  concern  that  different  cations  could  accumulate  differ¬ 
ent  amounts  of  hydration  water,  which,  in  turn,  affects  the  dielectric 
properties  of  the  sample,  homoionic  clays  could  be  used  in  future 
studies. 

d.  Soil  mixtures  containing  varying  amounts  of  a  pure  sand,  a  pure  silt, 
and  a  pure  clay  could  be  developed  and  tested  to  more  clearly  deter¬ 
mine  the  interrelationships  of  idealized  mixture  components  on  the 
dielectric  properties  of  the  mixture. 

e.  To  better  understand  the  effects  of  minerology,  monomineralic  samples 
could  be  prepared  for  testing. 

/.  All  of  the  measurements  made  in  this  study  were  done  on  samples  pre¬ 
pared  with  pure  distilled,  deionized  water.  The  effect  of  increased  ion 
availability  could  be  studied  by  preparing  samples  with  fresh  water  and 
with  bracWsh  water. 

g.  A  clearer  picture  of  the  low-frequency  loss  mechanisms  could  be  devel¬ 
oped  by  extending  these  measurements  to  lower  frequencies.  This 
would,  however,  require  a  new  set  of  instrumentation. 
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Appendix  A 

Data  Collection  and  Data 
Presentation  Software 


The  first  two  codes  contained  in  this  appendix  are  written  in  Hewlett 
Packard  BASIC  and  are  used  in  an  interactive  manner  on  an  HP9000/362  data 
controller.  The  first  code  instructs  an  HP8510C  Network  Analyzer  System  to 
collect  data,  stores  raw  data  (in  the  form  of  complex  S-parameters)  on  the 
HP9000/362,  and  produces  color,  annotated  graphs  of  the  measured  complex 
dielectric  constant  versus  frequency.  It  is  also  used  to  retrieve  existing  data 
files  for  further  study.  The  second  code  reads  existing  binary  S-parameter 
files  and  writes  ASCII  files  containing  real  and  imaginary  components  of  the 
complex  dielectric  constants  to  DOS-formatted  diskettes. 

The  other  codes  are  written  in  the  SAS  language  (SAS  Institute,  Inc. 

1988)^  and  are  used  for  creating  SAS  data  files  and  generating  the  plots  con¬ 
tained  in  this  report.  The  third  listing  is  the  program  that  creates  the  SAS 
data  set  that  provides  input  to  all  other  operations.  The  fourth  listing  is  the 
code  that  plots  data  and  calculated  results  versus  frequency  (Appendix  B),  and 
the  fifth  code  plots  data  and  calculated  results  versus  moisture  content 
(Appendix  C). 


'  References  cited  in  this  appendix  are  located  at  the  end  of  the  main  text. 
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Data  Collection  Program 


10  I  PROGRAM  NAME:  EPS12  3  MAR  94 
20  I 

30  I  A  MULTIPURPOSC  CODE  FOR  1)  COLLECTING  COAXIAL  REFLECTION  AND 
40  I  TRANSMISSION  bATA  AND  CALCULATING  BOTH  THE  REAL  AND  IMAGINARY 
50  I  PARTS  OF  THE  COMPLEX  DIELECTRIC  CONSTANT,  OR  2)  DISPLAYING 
60  !  THE  RESULTS  OF  PREVIOUSLY  COLLECTED  DATA 

70  ! 

80  I  USES  THE  HP9000  ACS  ALGORITHM  WITH  SIGN  ADJUSTMENTS 
90  I 

100  OPTION  BASE  1 

110  COMPLEX  S11{1:801),S21(l:801),Rhs,Lhs,Rt,Eps 

120  INTEGER  Preamble, Bytes,H!drnum,Temper,Hldr_tmp,Temper_tmp 

130  DIM  Ttl$I60J.Secllne$I801.Epsp(1:801),Epspp(r:801) 

140  ! 

150  GOSUB  Init^vars 
160  GOSUB  System^setup 
170  WHILE  1 
ISO  GOSUB  Main_menu 

190  GOSUB  Plot^data 

200  PAUSE 

210  LINPUT  "D  PRINTER  2)  PLOTTER  :  ",Plt$ 

220  IF  Plt$=’'2"  THEN 

230  GOSUB  Copy_data 

240  GOSUB  Plot^data 

250  END  IF 

260  IF  Plt$=:-r  THEN 

270  GOSUB  Plot^data 

280  GOSUB  Print_to_701 

290  GOSUB  Copy_data 

300  END  IF 

310  PRINTER  IS  CRT 

320  PLOTTER  IS  CRT.-INTERNAL" 

330  CLEAR  SCREEN 

340  KEY  LABELS  ON 

350  LOCAL  @Nwa 

360  END  WHILE 
370  STOP 

380  Msin_menu:  . ••*••••••*«••*•*•••• 

390  REPEAT 
400  PRINTER  IS  CRT 
410  BEEP2000,.2 
420  CLEAR  SCREEN 
430  PRINT 
440  PRINT 

450  PRINT  "  FILENAME  :",Namef$ 

460  PRINT  " _ 

470  PRINT 
480  PRINT 

490  PRINT**  N,TAU  =  ",Nlst.Tau 
500  PRINT 
510  PRINT 

520  PRINT  *•  MAX  &  MIN  VALUE  OF  EPSP  =  ".MAXIEpspC )).MIN(EpspC)) 

530  PRINT  "  MAX  &  MIN  VALUE  OF  EPSPP  =  ".MAX{Epspp(M).MIN(Epspp{M) 

540  PRINT 
550  PRINT 
560  PRINT 

570  PRINT  "  1)  ENTER  SAMPLE  HEADER  DATA  :-,Ttl$ 

580  PRINT  "  2)  CHANGE  PLOT  SCALE  -  MIN  MAX  TICK  :".Wmln.Wmax,Ytick 
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590  PRINT  -  3)  CHANGE  ELECTRICAL  DELAYS  -  S11S21  :",Edl  I^Edai/ps* 

600  PRINT  "  4)  CHANGE  CUT-OFF  VALUE  :",lmrhs 

610  PRINT  "  5)  CHANGE  SAMPLE  LENGTH 

620  PRINT  "  6)  FILE  MANAGEMENT  MENU  -  PERM  RTRV  :",Dirperm$,Dirrecv$ 
630  PRINT  “  7)  MAKE  ANOTHER  MEASUREMENT" 

640  PRINT  "  8)  EXIT  PROGRAM" 

650  Sell  =:0 

660  INPUT  "SELECT  :  ".Sell 
670  SELECT  Sell 
680  CASE  1 
690  GOSUBMenu3  * 

700  CASE  2 

710  INPUT  "INPUT  PLOT  SCALE  MINIMUM,  MAXIMUM,  AND  TICK  MARK  INCRE¬ 
MENTS "  ,Wmin.Wmax,Ytick 
720  CASE  3 

730  INPUT  "INPUT  ELECTRICAL  DELAYS  S11,S21  (ps)  :",Edn,Ed21 
740  OUTPUT  @Nwa;"Sl  l;ELED";Edl  1;"ps;" 

750  OUTPUT  @Nwa,*"S21;ELED";Ed21;"ps;" 

760  GOSUB  Get^data 

770  GOSUB  Compute  eps 

780  CASE  4 

790  INPUT  "INPUT  THE  CUTOFF  VALUE  :",lmrhs 
800  GOSUB  Get^data 

810  GOSUB  Compute_eps 

820  CASE  5 

830  INPUT  "ENTER  THE  SAMPLE  LENGTH  (cm)  7  ".D 

840  GOSUB  Get_data 

850  GOSUB  Compute_eps 

860  CASE  6 

870  GOSUB  Menu2 

880  CASE  7 

890  GOSUB  Get^newfile 

900  GOSUB  Get_data 

910  GOSUB  Compute_eps 

920  CASE  8 

930  STOP 

940  END  SELECT 

950  UNTIL  Sell  =0 

960  RETURN 

970  Menu2: 

980  REPEAT 
990  CLEAR  SCREEN 
1000  PRINT 

1010  PRINT  "  1)  RETRIEVE  AN  EXISTING  FILE 

1020  PRINT  "  2)  CHANGE  DIRECTORY  NAME  FOR  EXISTING  FILE  :",Dirrecv$ 

1030  PRINT 

1040  PRINT  "  3)  SAVE  FILE  IN  PERMANENT  DIRECTORY  :" 

1050  PRINT  "  4)  CHANGE  PERMANENT  DIRECTORY  NAME  :",Dirperm$ 

1060  PRINT  "  5)  CREATE  PERMANENT  DIRECTORY" 

1070  PRINT  "  6)  DELETE  FILE  FROM  PERMANENT  DIRECTORY" 

1080  PRINT 

1090  PRINT  "  7)  CLEAR  TEMPORARY  DIRECTORY  :  DATATEMP" 

1100  Sel2=0 

1110  INPUT  "SELECT  :  ",Sel2 
1120  SELECT  Sel2 
1130  CASE  1 
1140  GOSUB  Get_oldfile 
1150  GOSUB  Gerdata 
1160  GOSUB  Compute_eps 
1170  Sel2  =  0 
1180  CASE  2 
1190  CATDirroot$ 

1200  LINPUT  "ENTER  DIRECTORY  NAME  FOR  EXISTING  FILE  :".Dirrecv$ 

1210  CASE  3 
1  220  GOSUB  Copy  data 
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1 230  CASE  4 
1 240  CAT  Dirroot$ 

1 250  LINPUT  -ENTER  PERMANENT  DIRECTORY  NAME  :-,Dirperm$ 

1260  CASE  5 
1 270  CAT  Dirroot$ 

1 280  UNPUT  -ENTER  PERMANENT  DIRECTORY  NAME  ;-,Dirperm$ 

1 290  ON  ERROR  GOTO  1 270 
1 300  CREATE  DIR  Dirroot$&-r&Dirperm$ 

1310  CREATE  DIR  Dirroot$&-/-&Dirperm$&-/RAW" 

1320  OFF  ERROR 
1330  CASE  6 

1340  CAT  Dirroot$&-r&Dirperm$ 

1350  LINPUT  -ENTER  FILENAME  TO  DELETE  :-.Del$ 

1360  ON  ERROR  GOTO  1340 

1370  PURGE  Dirroot$&-r&Dirperm$&-r&Del$ 

1380  PURGE  Dirroot$&-/-&Dirperm$&-/RAWr&Del$&-R" 

1390  OFF  ERROR 
1400  CASE  7 

1410  CAT  Dirroot$&-/"&Dirtemp$ 

1420  LINPUT  -CLEAR  THIS  DIRECTORY  (Y/N)  :  -,Del$ 

1430  IF  Del$*-Y-  THEN 

1440  PURGE  Dirroot$&-/"&Dirtemp$&-/RAW/*" 

1450  PURGE  Dirroot$ &-r&Dirtemp $&"/•" 

1460  CREATE  DIR  Dirroot$&"r&Dirtemp$&"/RAW" 

1470  END  IF 
1480  END  SELECT 
1490  UNTIL  SeI2=0 
1500  RETURN 

1510  Menu3:  !•****•  . . . . . 

1520  REPEAT 
1 530  CLEAR  SCREEN 
1540  PRINT 

1550  PRINT  -  HEADER  DATA  :-,Ttl$ 

1560  PRINT _ _ _ "  ' 

1570  PRINT 

1580  PRINT- 1)  SAMPLE  :-.Samp$ 

1590  PRINT  -  2)  HOLDER  NUMBER  :-,Hldrnum 

1600  PRINT  -  3)  TEMPERATURE  (DEGREES  C)  Temper 
1610  PRINT 

1 620  PRINT  -  4)  INPUT  OWN  TITLE  - 
1630  Se(3=0 

1640  INPUT -SELECT  :  ",Sel3 
1650  SELECT  SeI3 
1660  CASE  1 

1 670  LINPUT  -INPUT  SAMPLE  LABEL  :",Samp$ 

1680  CASE  2 

1690  INPUT  -INPUT  HOLDER  NUMBER  :-,Hldrnum 
1700  CASE  3 

1710  INPUT  "INPUT  TEMPERATURE  (Degrees  C)  :",Temper 
1720  CASE  4 

1730  LINPUT  -INPUT  A  TITLE  FOR  THE  PLOT  :-.Ttl$ 

1740  Sel3=0 
1750  END  SELECT 

1760  IF  Sel3<>0  THEN  Ttl$  =  "Sample  "&Samp$&-  -  Holder#"&VAL$(Hldrnum)&-  * 
-&VAL${Temper)&-  DegC" 

1770  UNTIL  Sel3=0 
1780  RETURN 

1790  Get^newfile:  j« . . . . . 

1800  I 

1810  I  PREPARE  TO  MAKE  A  NEW  SET  OF  MEASUREMENTS 

1820  I 

1830  REPEAT 

1 840  Dup_filename  =*0 

1 850  Fd$  =  DATE$(TIMEDATE) 

1860  Ft$  -T1ME${TIMEDATE) 

1 870  Nameold$  =  Namef$ 
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1880  Namet$=Fd$11,21&UPC$(Fd$14,41)&UPC$(Fd$[6,61)&Fd$(11,111S<Ft$11,21&Ft$[4,51 

1890  Namer$  =Namef$&"R" 

1900  IF  Nameold$=Namef$  THEN 
1910  PRINT 

1920  PRINT  "DUPLICATE  FILES  EXIST  -  PLEASE  WAIT  !!!" 

1930  WAITS 
1 940  Dup^filename  =  1 

1950  END  IF 

1960  UNTIL  Dup_fiIename=0 
1970  CLEAR  SCREEN 
1980  PRINT 

1990  PRINT . MAKING  8510  MEASUREMENT . 

2000  PRINT 

2010  PRINT"  PLEASE  WAIT  mil" 

2020  OUTPUT  @Nwa;"WAIT;SING;" 

2030  OUTPUT  (9)Nwa;"WAIT;S11;ELED";Edl1;"ps;" 

2040  OUTPUT  @Nwa;"WAIT;S21  ;ELED";Ed21  :"ps;" 

2050  OUTPUT  @Nwa;"WAIT;Sl  1;F0RM3;0UTPDATA;" 

2060  ENTER  @Nwa_data;Preamble;Bytes 
2070  M  =  Bytes/16  “ 

2080  REDIM  Slid :M),S21  (1  :M) 

2090  ENTER  @Nwa_data;Sl  1  (*) 

2100  OUTPUT  @Nwa;"WAIT;S21;F0RM3;0UTPDATA;" 

2110  ENTER  @Nwa_data;Preamble;Bytes;S21  (•) 

2120  CREATE  BOAT  Dirroot$&"r&Dirtemp$&"r&Namef$d,HdrJgth  +  (Bytes"2) 

2130  ASSIGN  @F1  TO  Dirroot$&"r&Dirtemp$8i"r&Namef$;FORMAT  OFF 
2140  OUTPUT  @F1  USING  "10A";Namef$,Samp$ 

2150  OUTPUT  @F1  ;D,HIdrnunn, Temper, Edi  1  ,Ed21  ,lmrhs, Preamble, Bytes 
2160  OUTPUT  @Fl;Sll('').S21d) 

2170  ASSIGN  @F1  TO  • 

2180  OUTPUT  @Nwa;"WAlT;Sll;F0RM3;0UTPRAWl;" 

2190  ENTER  @Nwa_data;Preamble;Bytes,Sl  1  ("I 
2200  OUTPUT  @Nwa;"WAIT;S21  ;FORM3;OUTPRAW2;" 

2210  ENTER  @Nwa_data;Preamble;Bvtes;S2ir) 

2220  CREATE  BOAT  Dirroot$&"r&Dirtemp$&"/RAWr&Namer$.1  ,HdrJgth  + (Bytes*2) 
2230  ASSIGN  @F1  TO  Dirroot$&"r&Dirtemp$&"/RAWr&Namer$;FORMAT  OFF 
2240  OUTPUT  @F1  USING  "10A";Namer$,Samp$ 

2250  OUTPUT  @F1  ;D,Hldrnum, Temper, Edl  1  ,Ed21  .Imrhs, Preamble, Bytes 
2260  OUTPUT  @F1;S1 1{*).S21  (*) 

2270  ASSIGN  @F1  TO  " 

2280  RETURN 

2290  Get  oldfile:  ! . . 


!  RETRIEVE  AN  EXISTING  DATA  FILE  FOR  FURTHER  ANALYSIS 
I 

CAT  Dirroot$&"/"&Dirrecv$ 

BEEP  4000, .2 

INPUT  "WHAT  IS  THE  NAME  OF  THE  EXISTING  DATA  FILE?",Namef$ 

ON  ERROR  GOTO  2330 

ASSIGN  @F1  TO  Dirroot$&"/"&Dirrecv$&"r&Namef$;FORMAT  OFF 
OFF  ERROR 

ENTER  @F1  USING  -10A";Nameold$,Samp$ 

ENTER  @F1  ;D,Hldrnum,Temper,Ed1 1  ,Ed21 , Imrhs, Preamble, Bytes 

Ttl$  =  "Sample  "&Samp$&"  -  Holder#" &VAL$(Hldrnum)&"  -  "&VAL$(Temper)&" 

M  =  Bytes/1 6 
PRINT 

PRINT  "NUMBER  OF  DATA  POINTS  IN  FILE  =  ",M 
REDIM  S11{l:M),S21d:M) 

ENTER  @Fl;S11d),S21d) 

ASSIGN  @F1  TO  • 

OUTPUT  @Nwa;"WAIT;S1 1  ;ELED";Ed1 1  ;"ps;" 

OUTPUT  @Nwa;"WAlT;S21;ELED";Ed21;"ps;" 

Last  =  LEN(Namef$) 

IF  Namef$lLast,Lastl  =  "R"  THEN 


Appendix  A  Data  Collection/Data  Presentation  Software 


2520  Namer$  »Namef$ 

2530  Namef$lLast,LastI  "" 

2540  Namef$  *TRlM$(Namef$) 

2550  PRINT 

2560  PRINT  "  LOADING  RAW  DATA  RLE:  ",Namer$ 

2570  PRINT  "MUST  LOAD  CORRESPONDING  CALIBRATION  RLE  ON  8510" 

2580  PRINT 

2590  PRINT  "  —  PRESS  CONTINUE  WHEN  READY - " 

2600  PAUSE 

2610  PRINT 

2620  PRINT  "LOADING  RAW  DATA  RLE  . " 

2630  OUTPUT  @Nwa;"WAIT;HOLD;S1 1;FORM3;INPURAW1;" 

2640  OUTPUT  @Nwa_data;Pr0annbIe, Bytes, S1 1  (*) 

2650  OUTPUT  (5)Nwa;"WAIT;S21  ;F0RM3;INPURAW2;" 

2660  OUTPUT  @Nwa_data;Preamble, Bytes, S21  (*) 

2670  ELSE 

2680  OUTPUT  @Nwa;"WAIT;HOLD;Sl  1  ;FORM3;INPUDATA;" 

2690  OUTPUT  @Nwa__data;Preamble,Bytes,S1 1  (*) 

2700  OUTPUT  @Nwa;"WAIT;S21  ;F0RM3;INPUDATA;" 

2710  OUTPUT  @Nwa_data;Preamble,Bytes,S21  (*) 

2720  END  IF 

2730  ASSIGN  @F1  TO  • 

2740  RETURN 

2750  Get_data:  (*•••*••••••*••*'•« . . 

2760  I  ESTIMATE  THE  INITIAL  'N'  VALUE  FOR  COMPUTING  THE  CORRECT  PHASE  ANGLE 
2770  PLOTTER  IS  CRT,"1NTERNAL" 

2780  OUTPUT  @Nwa;"WAIT;Sl  1  ;WAIT;AUTO;" 

2790  OUTPUT  @Nwa;"S21;DELA;" 

2800  OUTPUT  @Nwa;"WAIT;MARKr;Fi;";" 

2810  OUTPUT  @Nwa;"OUTPMARK;" 

2820  ENTER  @Nwa;Tau, Dummy 
2830  N  =  INT(Tau*Fi) 

2840  IF  N<0THEN  N=0 
2850  IF  N>99  THEN  N=99 

2860  OUTPUT  @Nwa;"WAlT;REIP;WAIT;AUTO;ENTO" 

2870  I 

2880  I  LOAD  THE  DATA  INTO  COMPLEX  ARRAYS  ON  THE  330  FOR  FURTHER  PROCESSING 
2890  I 

2900  OUTPUT  @Nwa;"WAIT;S11;FORM3;OUTPFORM;" 

2910  ENTER  @Nwa_data;PreambIe;Bytes:S1 1  (*) 

2920  OUTPUT  @Nwa;"WAIT;S21  ;FORM3;OUTPFORM;" 

2930  ENTER  @Nwa_data;Preamble;Bytes;S21  (*) 

2940  RETURN 

2950  Compute_eps:  i*»***»««»«* . . . 

2960  I  CALCULATE  THE  COMPLEX  DIELECTRIC  CONSTANT 
2970  C  =  3.E  +  10 
2980  RAD 

2990  PRINTER  IS  CRT 
3000  N1st  =  N 
3010  Anglel=0. 

3020  FOR  1  =  1  TO  Np 

3030  Rhs  =  (1  +S21(l)‘‘S21(!)-Sn(l)-Sn(l))/(2*S21(!)) 

3040  IF  ABS(Rhs)>1.E  +  6  THEN  GOTO  3220 

3050  IF  ABS(IMAG(Rhs))<lmrhsTHEN  Rhs  =  CMPLX(REAL(Rhs),0) 

3060  Lhs  =  ACS(Rhs) 

3070  Lhs  =  CMPLX{REAL(Lhs).ABS(lMAG(Lhs))) 

3080  IF  IMAG(Rhs)<0  THEN  Lhs  =CMPLX(2‘Pl-REAL(Lhs).!MAG(Lhs)) 

3090  Angle  =  REAL(Lhs) 

3100  IF  1  =  1  THEN  GOTO  3130 

3110  IF  ABS{Angle-Anglel)<4.5  THEN  GOTO  3130 

3120  N  =  N  +  1 

3130  Lhs  =  CMPLX(REAL(Lhs)  +  N*2"PUMAG(Lhs)} 

3140  I  PRINT  USING  "DDDD.5{2X.SDDD.DDD). 2X,0DD";I.Rhs,Lhs, Angle, N 
3150  Anglel  =Angle 
3160  Rhs  =  Lhs"Lhs 

3 1 70  Radf  req  =  2  *  PI  *  (Fi  +  (I- 1 }  *  (Ff-Fi)/(Np- 1 )) 
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3180  Fact«C/(Radfreq*D) 

3190  Epsp(l)- Fact  “Fact*  REAL(Rhs) 

3200  Epspp(l)  =  Fact»Fact*IMAG(Rhs) 

3210  I  PRINT  Epsp(l),Epspp(l) 

3220  NEXTl 
3230  RETURN 

3240  Plot^data:  I*** . . . 

3250  I 

3260  !  PLOT  THE  RESULTS 

3270  I 

3280  GINIT 

3290  PRINTER  IS  CRT 

3300  PLOTTER  IS  CRT, "INTERNAL" 

3310  IFPIt$="2"THEN 
3320  PRINTER  IS  705 
3330  PLOTTER  IS  705,"HPGL" 

3340  END  IF 
3350  GRAPHICS  ON 
3360  CLEAR  SCREEN 
3370  Plt$="0" 

3380  X_gdu_max  =  100»MAX(1, RATIO) 

3390  Y~gdu_max  =  100*MAX(1,1/RAT10) 

3 400  Xmax  =  X_gdu_max  *1.0 
3410  Ymax=Y_gdu_max*1 .0 
3420  LORG  6 
3430  CSI2E3.0 
3440  PEN  6 

3450  MOVE  Xmax*.5,Ymax 
3460  LABEL  Ttl$ 

3470  Secline$="FILE:  "&Namef$&"  SAMPLGTH  =  "&VAL$(D)&"cm  DE- 
LA  YS(S1 1  ,S21 )  =  "&VAL$(Ed1 1  )&"."&VAL${Ed21  )&"ps" 

3480  MOVE  Xmax*.5,YmaxV97 
3490  LABEL  Secline$ 

3500  Secllne$="lmrhs  =  "&VAL$(lmrhs)&"  N  =  "&VAL${N1st)&"  Tau  =  "&VAL$(Tau/1.0E-12)- 
&"ps  "&Dirperm$ 

3510  MOVE  Xmax*.5.Ymax*.94 
3520  LABEL  Seciine$ 

3525  MOVE  XmaxV5..51  *Ymax 

3526  DIMSpar$(80] 

3527  Spar$="S11  at  45  MH2=  "&VAL$(REAL(S1 1  (1)))&"  .  "&VAL$(1MAG(S1 1  (1))) 

3528  LABEL  Spar$ 

3529  MOVE  XmaxV5,.48*Ymax 

3531  Spar$="S21  at45MH2=  "&VAL$(REAL{S21  {1)))&"  .  "&VAL$(IMAG(S21  (1))) 

3532  LABEL  Spar$ 

3533  DEG 
3540  LDIR  90 

3550  MOVE  Xmax^04.Ymax^75 
3560  LABEL  "EPS'" 

3570  MOVE  Xmax •.04, YmaxV3 
3580  LABEL  "EPS"" 

3590  LORG  4 

3600  LDIRO 

3610  KEY  LABELS  OFF 

3620  MOVE  Xmax *.5,0 

3630  LABEL  "FREQUENCY  (GH2)" 

3640  VIEWPORT  .15*Xmax,.9*Xmax..55*Ymax..90*Ymax 
3650  GOSUB  Plot_axes 
3660  FOR  1  =  1  TO  Np 

3670  IF  ABS(Epsp{l))<.001  THEN  GOTO  3710 
3680  X  =  LGT(Fi  +  {I- 1 )  *  (Ff-Fi)/(Np- 1 )) 

3690  PEN  2 
3700  PLOT  X,Epsp(l) 

3710  NEXTl 

3770  VIEWPORT  .1 5* Xmax,. 9 "Xmax,.  1  •  Ymax,.45« Ymax 
3780  GOSUB  Plot  axes 
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3790  FOR  1*1  TO  Np 

3800  X  *  LGT(R  +  (I-  D  MFf-Fi}/(Np- 1 )) 

3810  PEN  2 

3820  PLOT  X,Epspp(l) 

3830  NEXT  I 
3840  RETURN 

3850  Plot_ax8s:  ••••**••■  * . 

3860  WINDOW  7,1 1,Wmin,Wmax 
3870  PEN  5 

3880  AXES  1,Ytick,7,0r1, 1,2 
3890  UNE  TYPE  4.2  * 

3900  GRID  1,Ytick,7, 0,1, 1,2 

3910  LINE  TYPE  1.5 

3920  CLIP  OFF 

3930  CSIZE  2.5,.5 

3940  LORG  6 

3950  FOR  1  =  8  TO  11 

3960  MOVE  l.-.2*Ytick 

3970  LABEL  USING  "#,K";(10"I)/1  .E  +  9 

3980  F0RJ  =  2T0  9 

3990  Jtick  =  LGT(J) 

4000  MOVE  1-1  +Jtick,Wmin 
4010  PEN  5 

4020  PLOT  1-1  -I- Jtick,.04*(Wmax-Wmln),-1 
4030  NEXT  J 
4040  NEXT  I 
4050  LORG  8 

4060  FOR  I  =  Wmin  TO  Wmax  STEP  Ytick 

4070  MOVE  6.85,1 

4080  LABEL  USING  "#,SDDD.D";) 

4090  NEXT  I 
4100  PENUP 
4110  RETURN 

4120  Print_to_701:  I . * . . . . . . 

4130  I  THESE  4  LINES  MARKS  START  OF  NOISE  &  PRINTS  8510  SCREEN  TO  XL300. 
4140  OUTPUT  @NwarWAIT;OUTPERRO;" 

4150  ENTER  @Nwa_data;Preamble 

41 60  OUTPUT  @Nwa;"WAlT;S1 2;WA1T;AUT0;WA1T;MARK1  ;TARV  -80;MARKTARG;" 
4170  OUTPUT  (5)Nwa;"WAlT;ENT0;PRlNALL;" 

4180  LOADSUB  ALL  FROM  "/UTILITIES/GDUMP^C" 

4190  Gdump_colored(CRT.701."N0RMALM80,"OFF") 

4200  PRINTER  IS  701 
4210  PR1NTCHR$(12) 

4220  PRINTER  IS  CRT 
4230  Plt$="0" 

4240  BEEP  2000.. 2 

4250  PRINT  -PRINT  COMPLETE" 

4260  RETURN 

4270  CopY_data:  . . . . . . 

4280  OUTPUT  @Nwa;"WAIT;S1 1  ;F0RM3;0UTPDATA;- 
4290  ENTER  @Nwa_data;Preamble,Bytes,Sl  1  (*) 

4300  OUTPUT  (s)Nwa;"WAlT;S21;F0RM3;0UTPDATA;- 
4310  ENTER  @Nwa_data;Preamble,Bytes,S21  (•) 

4320  Namer$  =Namef$&"R" 

4330  ON  ERROR  GOTO  4360 

4340  PURGE  Dirroot$&"r&Dirperm$&"/-&Namef$ 

4350  PURGE  Dirroot$&"r&Dirperm$&"/RAWr&Namer$ 

4360  OFF  ERROR 

4370  CREATE  BOAT  Dirroot$&-/"&Dirperm$&-r&Namef$,1  .Hdrjgth-H  (Bytes*2) 

4380  ASSIGN  @F1  TO  Dirroot$&"r&Dirperm$&-r&Namef$;FORMAT  OFF 
4390  OUTPUT  @F1  USING  -10A-;Namef$.Samp$ 

4400  OUTPUT  @F1  ;D,Hldrnum,Temper,Edl  1  ,Ed21  .Imrhs.Preamble.Bytes 
4410  OUTPUT  @F1;S11(*).S21C) 

4420  OUTPUT  @Nwa;"WAIT;Sl  1  ;F0RM3;0UTPRAW1 
4430  ENTER  @Nwa  datajPreamble.Bytes.S1 1  (*) 

4440  OUTPUT  @Nwa;"WAIT;S21  ;FORM3;OUTPRAW2;- 
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4450  ENTER  @Nwa_data;Preamble,Bytes,S21  {*) 

4460  CREATE  BOAT  Dirroot$&"r&Dirperm$&"/RAWr&Namer$,l,HdrJoth  +  (Bytes*2) 
4470  ASSIGN  TO  Dirroot$&"r&Dirperm$&’*/RAW/"&Namer$;FORMAT  OFF 

4480  OUTPUT  @F1  USING  "10A'';Nanner$,Samp$ 

4490  OUTPUT  @F1  ;D,Hldrnum,Temper,Ed1 1  ,Ed21  .Imrhs.Preamble.Bytes 
4500  OUTPUT  @F1;S11C),S21C) 

4510  ASSIGN  @F1  TO  • 

4520  RETURN 
4530  System^setup: 

4540  WILDCARDS  DOS 

4550  REDIM  SI  1  (1  :Num_data_pts),S21  (1  :Num_data_pts) 

4560  FOR  I  =  1  TO  Num_data_pts 
4570  S21  (0  =  1.0 

4580  NEXTl 

4590  ASSIGN  @Nwa  TO  716 

4600  ASSIGN  @Nwa_dataTO  716;FORMAT  OFF 

4610  LOCAL  @Nwa 

4620  PRINTER  IS  CRT 

4630  GINIT 

4640  PLOTTER  IS  CRT/INTERNAL" 

4650  GRAPHICS  ON 
4660  ASSIGN  (s>Hpib  TO  7 
4670  REMOTE  (9>Hpib 
4680  ABORT  <s)Hpib 
4690  CLEAR  @Nwa 
4700  CLEAR  71 7 
4710  CLEAR  SCREEN 
4720  BEEP  2000, .2 

4730  PRINT  "  NUMBER  OF  DATA  POINTS  =",Num_data_,pts 

4740  INPUT  "ENTER  INSTRUMENT  STATE  NUMBER  TO  RECALL  :",lnstr_num 

4750  IF  lnstr_num<>0  THEN 

4760  PRINT 

4770  PRINT  "WAIT  -  MAKING  INITIAL  MEASUREMENT  FOR  INSTr.lnstr_num 
4780  OUTPUT  @Nwa;"WAIT;RECA"&VAL$((nstr_num)&";" 

4790  OUTPUT  <slNwa;"WAIT;Sl  1  ;ELED";Edl  l;"ps;" 

4800  OUTPUT  @Nwa;"WAIT;S21  ;ELED";Ed21  ;"ps;WAIT;ENTO;" 

4810  OUTPUT  @Nwa;"WAIT;SING;WAIT;HOLD;" 

4820  END  IF 
4830  RETURN 

4840  lnit_vars:  }••••*•**•••••**•*•••*•**•■•••*•••«*•••••••*•• 

4850  Num_data__pts  =  801  !  #  OF  8510  DATA  POINTS 

4860  PRINT  "NUMBER  OF  DATA  POINTS  IS",Num_data_pts 
4870  INPUT  "INPUT  A  NEW  NUMBER  OF  POINTS",Num_^data_pts 
4880  Samp$=" AIR  "  I  SAMPLE  NAME 

4890  PRINT  "MEASUREMENT  DESCRIPTION  lS",Samp$ 

4900  LINPUT  "INPUT  A  NEW  DESCRIPTION",Samp$ 

4910  D  =  10.0  I  HOLDER  LENGTH 

4920  PRINT  "HOLDER  LENGTH  IS",D."  CM" 

4930  INPUT  "INPUT  A  NEW  HOLDER  LENGTH". D 

4940  Hldrnum  =  2  !  HOLDER  NUMBER 

4950  PRINT  "HOLDER  NUMBER  IS".HIdrnum 

4960  INPUT  "INPUT  A  NEW  HOLDER  NUMBER", HIdrnum 

4970  Temper  =  22  !  TEMPERATURE 

4980  PRINT  "SAMPLE  TEMPERATURE  IS",Temper."  DEG  C" 

4990  INPUT  "INPUT  A  NEW  SAMPLE  TEMPERATURE", Temper 

5000  Edi  1=270.0  I  S1 1  DELAY 

5010  PRINT  "S1 1  DELAY  IS".Ed  11,"  PICOSECONDS" 

5020  INPUT  "INPUT  A  NEW  S1 1  DELAY", Edi  1 

5030  Ed21  =270.0  I  S21  DELAY 

5040  PRINT  "S21  DELAY  IS".Ed21,"  PICOSECONDS" 

5050  INPUT  "INPUT  A  NEW  S21  DELAY", Ed21 
5060  lmrhs=0.  !  CUT-OFF  VALUE 

5070  Ttl$  =  "Sample  "&Samp$&"  -  Holder#"&VAL$(Hldrnum)&"  -  "&VAL$(Temper)&" 

DegC"  !  TITLE  ON  PLOTS 

5080  Hdr  stringl  =10  !  EVEN  NUMBER  OF  CHARS 


Appendix  A  Data  Collection/Data  Presentation  Software 


5090  Hdr_strino2«10  I  EVEN  NUMBER  OF  CHARS 
5100  Hdr_strings  *=  (4  +  Hdr_stringl )  +  (4  +  Hdr_string2) 

51 10  Hdr  integers  =  4  I  #  OF  INTEGRES  IN  HEADER 

5120  HdrJloats  =  4  I  it  OF  FLOATS  IN  HEADER 

5130  Hdr  lgth  =  Hdr_strings  +  (HdrJntegers*2)  +  (Hdr_floats*8) 

5140  Wfnfn=0  I  MINIMUM  PLOT  SCALE 

5 1 50  Wmax  =  2.0  I  MAXIMUM  PLOT  SCALE 

5160  Ytick  =  .5  I  TICK  MARKS 

5170  Dirroot$  ="/EPS/DATA"  I  ROOT  DIRECTORY  OF  DATA 

5180  Dlrperm$  =-PLOT;S"  I  PERMANENT  DIRECTORY.SAVED  WHEN  PLOTTED 

51 90  Dirr®cv$  =-DATATEMP"  I  DIRECTORY  OF  EXISTING  FILES  TO  RECEIVE 

5200  Dirtemp$  =*DATATEMP"  I  ALL  MEASUREMENTS  MADE  SAVED  HERE 

5210  Fi«4.5E  +  7  I  START  FREQUENCY 

5220  Ff  =  2.65E  +  10  I  STOP  FREQUENCY 

5230  Np  s=Num_dataj3ts  I  U  OF  8510  DATA  POINTS 

5240  RETURN 

5250  Errs:  PRINT . * . * . 

5260  PRINT  "  UNEXPECTED  ERROR  #**,ERRN."  AT  LINE#".ERRLN 
5270  PRINT  •  PRESS  CONTINUE  OR  RESET" 

5280  PAUSE 
5290  ERROR  RETURN 
5300  END 
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Program  for  Producing  ASCII  Data  Files 


10  I  PROGRAM  MAKEASCII  -  THIS  PROGRAM  READS  S1 1  &  S21  CORRECTED 
20  !  DATA  FILES,  CALCULATES  EP'  &  EP'',  AND 

30  I  STORES  DATA  IN  ASCII  TO  DOS  FLOPPY  DISK. 

40  OPTION  BASE  1 

50  COMPLEX  SI  1  (1:801), S21(l:801).Rhs,Lhs,Rt,Eps 

60  INTEGER  Preamble,Bytes,Hldrnum,Tempor,Phscorr,Statl,Stat2 
70  REAL  Fq,Fq1  ,Fq2,Fqstep,Fq1  lg,Fq2lg,Steplg 

80  DIM  TtI$[601,Secline$I60],Epsp{l:801),Epspp{l:801).ldx(l:801).Ep(1:801),Epp{l:801),Fr- 

eq{l:801) 

90  I 

100  WILDCARDS  DOS 

1 1 0  Instr^num  =  1  I  INSTRUMENT  STATE  NUMBER 

120  NurrTdata^pts^SOI  I  NUMBER  OF  DATA  POINTS 
130  IF  Num_dMa__pts<>801  THEN 
140  REDIM  Slid :Num_data_pts).S21  (1  :Num__data_pts) 

1 50  REDIM  Epspd  :Num_data_pts),Epsppd  :Num_data_pts) 

160  END  IF 

170  D  =  10.0 

180  Edl1=270 

190  Ed21=270 

200  lmrhs=0 

210  Fi  =  4.5E  +  7 

220  Ff  =  2.65E+10 

230  Np  =  Num_data_pts 

240  ASSIGN  @Nwa  TO  71 6 

250  ASSIGN  @Nwa_data  TO  716;FORMAT  OFF 

260  LOCAL  (iJNwa 

270  PRINTER  IS  CRT 

280  GINIT 

290  PLOTTER  IS  CRT, "INTERNAL" 

300  GRAPHICS  ON 
310  TotaIj3ts  =  100 

320  REDIM  Idxd  :TotaI_pts),Ep(1  :Total_pts).Epp(l  :Total_pts),Freq(1  :Total_pts) 

330  Fql  =  1 .00399375E  +  9 

340  Fq2  =  2.65E  +  10 

350  Fqstep  =  (Ff-Fi)/(Num_data_pts-1 ) 

360  ldx_start  =  (Fql  >Fi)/Fqstep  +  1 
370  Ptslg=TotaI_pts-ldx_start 
380  Fq1Ig  =  LGT(Fq1) 

390  Fq2lg  =  LGT(Fq2) 

400  Steplg  =  (Fq2lg-Fq1lg)/Ptslg 

41 0  FOR  I  =  1  TO  ldx_start 

420  Idxd)  =  I 

430  NEXT  I 

440  J  =  ldx_start 

450  K  =  ldx_start  +  1 

460  FOR  Fq  =  Fq1  TO  Fq2  STEP  Fqstep 

470  J=J  +  1 

480  IF  (10^(Fqllg  +  (Steplg*(K-ldx_start)))<Fq)  THEN 

490  ldx(K)=J*l 

500  K  =  K+1 

510  END  IF 

520  NEXTFq 

530  I 

540  ASSIGN  @Hp\b  TO  7 
550  REMOTE  @Hp\b 
560  ABORT  @Hpib 
570  CLEAR  (3)Nwa 
580  CLEAR  717 
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590  OUTPUT  @Nwa;"WAIT;RECA-&VAL$(lnstr^num)&";"  \  RECALL  INSTRUMENT  STATE 

600  I - 

610  CLEAR  SCREEN 

620  UNPUT  "ENTER  FILENAME  FOR  DATA  TO  DISKETTE  :  ",Ascname$ 

621  LINPUT  "ENTER  FILENAME  FOR  FREQUENCIES  TO  DISKETTE  :",Frqname$ 

630  CREATE  Ascname$&":DOS,1400",l 

640  ASSIGN  @F2  TO  Ascname$&":DOS,1400";FORMAT  OFF 

641  CREATE  Frqname$&":DOS,1400",l 

642  ASSIGN  @F3  TO  Frqname$&":DOS.1400";FORMAT  OFF 
650  Cntsaved^^O  ^ 

660  I - - - 

670  BEEP4000,.2 
680  CLEAR  SCREEN 
690  CAT  "/EPS/DATA" 

700  LINPUT  "ENTER  DIRECTORY  NAME  FROM  ABOVE  :  ",Dirname$ 

710  CLEAR  SCREEN 
720  ON  ERROR  GOTO  1 470 

730  CAT  "/EPS/DATAr&Dirname$;NO  HEADER.COUNT  Numjiles 
740  OFF  ERROR 

770  ALLOCATE  Filebuffer${l  :Num_files)I20] 

780  CAT  "/EPS/DATA/"&Dirname$  TO  Filebuffer$(*);NAMES 

781  MAT  SORT  Filebuffer$C} 

790  CLEAR  SCREEN 

791  PRINT  USING  "20A,/";Filebuffer$C) 

796  PRINT  "PRESS  CONTINUE" 

797  PAUSE 

800  FOR  1  =  1  TO  2 

810  IF  1  =  2  THEN  PRINTER  IS  701 

820  PRINT  USING 

830  PRINT  USING  "14A.r;Ascname$ 

840  PRINT  "DIR  FILE  SAMP  LGTH  #  MC  TMP" 

850  PRINT  " - - 

860  PRINTER  IS  CRT 
870  NEXT  I 

880  FOR  Cm  =  1  TO  Numjiles 
890  BEEP  2000,. 2 

900  Fmtl:  IMAGE  2D,2X.4A,2D,4X,1 2A.14A.# 

910  PRINTER  IS  701 

920  PRINT  USING  Fmtl  ;Cnt,"OF",Num  files, Dirname$,Filebuffer$(Cnt) 

930  PRINTER  IS  CRT 

940  PRINT  USING  Fmtl :Cnt, "OF", Numjiles, Dirname$,FiIebuffer$(Cnt) 

941  IF  Filebuffer$(Cnt)ll,31  =  "RAW"  THEN  GOTO  2100 

950  I  LINPUT  "CONVERT  FILE  (N=NO  Q  =  QUIT  B  =  BACK)  ?  ",Scl$ 

960  I  IF  Sel$  ="Q-  THEN  STOP 

970  I  IF  Sel$  =  "N"  OR  Sel$="B"  THEN 

980  I  PRINTER  IS  701 

990  I  PRINT 

1000  !  PRINTER  IS  CRT 

1010  I  PRINT 

1020  I  IF  SeI$="B"  THEN  Cnt  =  Cnt-2 
1030  I  GOTO  2100 
1040  !  END  IF 

1 050  Cntsaved  =  Cntsaved  +  1 
1060  Namef$  =sFilebuffer${Cnt) 

1070  ON  ERROR  GOTO  Errs 

1080  ASSIGN  @F1  TO  "/EPS/DATA/"&Dirname$&"r&Namef$;FORMAT  OFF 

1090  OFF  ERROR 

1091  Last_samp$  =  Samp$ 

1100  ENTER  @F1  USING  "10A";Namef$,Samp$ 

1110  ENTER  @F1  ;D,Hldrnum,Temper,Ed1 1  ,Ed21  .Imrhs, Preamble, Bytes 

1120  M  =  Bytes/16 

1130  REDIM  S11(1:M),S21(1:M) 

1140  ENTER  @F1;S11 0.521  (•) 

1150  ASSIGN  @F1  TO* 

1 1 60  Fmt4:  IMAGE  8A.2D.b,# 

1170  PRINTER  IS  701 
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1180  PRINT  USING  Fmt4;Samp$,D 
1190  PRINTER  IS  CRT 

1200  PRINT  USING  Fmt4;Samp$,D 

1201  BEEP  3000,. 2 

1202  INPUT  "MC?",Mc 

1210  !  IF  D  =  5.0  OR  Samp$(1 ,31  =  "AIR"  THEN 

1220  I  Mc=0 

1221  I  Last_sannp_air  ==  1 
1230  lELSE 

1240  I  IF  Samp$|l,5I»:"WATER*'  OR  Samp$[l,31  =  "H20-  THEN 
1250  I  Me  =  100 

1260  lELSE 

1261  !  IF  Last_samp_air  OR  Last_samp$  <  >  Samp$  THEN 
1263  I  BEEP3000..2 

1300  I  .  INPUT  "MC  FOR  1.4cm?",Mc1 
1303  I  IF  Mel  =999  THEN  STOP 
1310  I  INPUT  "MC  FOR  9.7em?",Mc2 
1341  I  END  IF 

1343  I  IF  D  =  1.4  OR  D  =  1.5  THEN  Me  =  Mel 
1345  I  IF  D  =  9.7  OR  D  =  10.0  THEN  Mc  =  Me2 

1350  1  END  IF 

1351  I  Last_samp_air  =  0 
1360  I  END  IF 

1370  Fmt2:  IMAGE  6D,7D.D,7D.4D 
1380  PRINTER  IS  701 

1390  PRINT  USING  Fmt2;Hldrnum,Mc,Temper,Cntsaved 
1400  PRINTER  IS  CRT 

1410  PRINT  USING  Fmt2;Hldrnum,Mc,Temper,Cntsaved 
1420  OUTPUT  @Nwa;"WAIT;H0LD;S1 1;F0RM3;INPUDATA;" 

1430  OUTPUT  @Nwa_data.*Preamble,Bytes.S1  KM 
1440  OUTPUT  @Nwa?’WAlT;S21;FORM3;lNPUDATA;" 

1450  OUTPUT  @Nwa_data;Preamble,Bytes,S2K*) 

1460  I 

1470  I  ESTIMATE  THE  INITIAL  'N'  VALUE  FOR  COMPUTING  THE  CORRECT  PHASE  ANGLE 
1480  1 

1490  PLOTTER  IS  CRT/INTERNAL** 

1500  OUTPUT  @Nwa;" WAIT;S  11  ;WAIT; AUTO;" 

1510  OUTPUT  @Nwa;"S21;DELA;" 

1520  OUTPUT  @Nwa;"WAIT;MARK1";Fi;";" 

1530  OUTPUT  @Nwa;"OUTPMARK;" 

1540  ENTER  @Nwa;Tau, Dummy 
1550  N  =  INT(Tau*Fi) 

1560  IF  N<0THEN  N=0 
1570  IF  N>99  THEN  N  =  99 

1580  OUTPUT  @>Nwa;"WAIT;REIP;WAIT;AUTO;ENTO" 

1590  I 

1600  I  LOAD  THE  DATA  INTO  COMPLEX  ARRAYS  ON  THE  330  FOR  FURTHER  PROCESSING 
1610  I 

1620  OUTPUT  @Nwa;"WAIT;S  11  ;FORM3;OUTPFORM;" 

1630  ENTER  @Nwa_data;PreambIe;Bytes;S1  K*) 

1 640  OUTPUT  @Nwa;"WAIT;S21  ;FORM3;OUTPFORM;" 

1650  ENTER  @Nwa_data;Preamble;Bytes;S2K*) 

1660  I 

1670  I  CALCULATE  THE  COMPLEX  DIELECTRIC  CONSTANT 
1680  1 

1690  C  =  3.E  +  10 
1700  RAD 

1710  PRINTER  IS  CRT 
1720  Nlst  =  N 

1730  Anglel=0. 

1731  PRINT  "WARNING  I!!!!!!!!!!!!!!" 

1732  PRINT  "RESETTING  SAMPLE  LENGTHS  TO  1 .4  cm  AND  9.7  cm" 

1733  IF  D  =  1.5  THEN  D  =  1.4 

1734  IF  D  =  10.  THEN  D  =  9.7 
1740  FOR  1  =  1  TO  Np 
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1750  Rhs  =  {l  +S21{l)»S21(l)-Sn(l)*Sll(l))/(2‘'S2im) 

1760  IF  ABS(Rhs)>1.E  +  6THEN  GOTO  1940 

1770  IF  ABS(IMAG{Rhs))<lmrhs  THEN  Rhs  =  CMPLX{REAURhs),0) 

1780  Lhs«ACS(Rhs) 

1790  Lhs=CMPLX(REAL(Lhs).ABS{IMAG{Lhs))) 

1800  IF  IMAGIRhsXO  THEN  Lhs  =  CMPLX{2*Pi-REAULhs),IMAG{Lhs)) 

1810  Angles  RE  AULhs) 

1 820  IF  I  =  1  THEN  GOTO  1 850 

1830  IF  ABS(Angle-Anglel)< 4.5  THEN  GOTO  1850 

1840  N=N  +  1  • 

1850  Lhs=CMPLX{REAL(Lhs)  +  N*2*PMMAG{Lhs)) 

1860  I  PRINT  USING  •DDDD,5(2X,SDDD.D0D),2X.DDD";l,Rhs,Lhs.Angle,N 
1 870  Angle!  =  Angle 
1880  Rhs  =  Lhs*Lhs 

1 890  Radf req  =  2  •  PI  *  {Fi  +  (1-1 )  •  (Ff-Fi )/(Np- 1 )) 

1 900  Fact  =  C/(Radf  req  •  D) 

1910  Epsp(l)  =  Fact*Fact*REAL(Rhs) 

1920  Epspp(l)  sFact*Fact*IMAG{Rhs) 

1930  I  PRINT  l,Epsp(l).Epspp(l} 

1940  NEXT  I 
1950  i 

1 960  FOR  I  =  1  TO  Total__pts 
1 962  Freq(l)  =  (Fi  +  ((ldx(l)-l)*Fqstep))/l  .OE  +  9 
1970  Ep(!)=Epsp(ldx(l)) 

1980  Epp(l)sEpspp{ldx(l)) 

1990  IF  1  =  1  THEN  GOTO  2020 

1991  IF  Ep{l)>999  THEN  Ep(l)=999 

1992  IF  Epp(l)  > 999  THEN  Eppd)  =  999 

2000  I  IF  ABS(Ep(l)-Ep(I-l))>2  THEN  Ep(l}  =  999 
2010  I  IF  ABS(Epp(l)-Epp(l-1))>2  THEN  Eppd)  =999 
2020  I 

2030  I  PRINTER  IS  701 
2040  !  PRINT  ldxd).I.Epd).Eppd) 

2050  I  PRINTER  IS  CRT 
2060  NEXT  I 

2070  Fmt3:  IMAGE  K./.K,/,K./,K,/.K,/,K 

2080  OUTPUT  @F2  USING  Fmt3;Namef$.Samp$,D,Hldrnunn,Mc, Temper 
2090  OUTPUT  @F2  USING  "3D.DD,r;EpC).Epp{*) 

2095  IF  Cntsaved  =  1  THEN  OUTPUT  USING  "2D.3D.r;Freq{*) 

2100  NEXTCnt 
2110  PRINT 

2120  PRINT  "FINISHED  DIRECTORY  ",Dirname$,"  PRESS  CONTINUE  FOR  MORE  DIRECTO¬ 
RIES" 

2130  PAUSE 

2140  DEALLOCATE  Filebuffer$C) 

2150  GOTO  670 
2160  STOP 
2170  Errs:  I 

2180  PRINT . . . . . . 

2190  IF  ERRN  =  54THEN 

2200  PRINT  "  NO  COPY  -  FILE  ALREADY  EXIST!" 

2210  ELSE 

2220  PRINT  "  UNEXPECTED  ERROR  #".ERRN."  AT  LINE#",ERRLN 
2230  END  IF 

2240  PRINT  "  PRESS  CONTINUE  OR  RESET" 

2250  PAUSE 
2260  ERROR  RETURN 
2270  END 
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SAS  Data  File  Creation  Code 


*  program  'makessdB* 

• . 

lib  name  eps  'c;\sas_win\aug93\'; 
data  eps.aug93_5; 
infiie  '*aug93\aug93_5^ 

input  date  $10.  material  $10.  holder  hidrnum  moisture 
temp  el -el  00  epI-eplOO; 
array  e{100)  el -el  00; 
array  ep{100)  epI-eplOO; 
do  i  =  1  to  1 00; 
if  e{i)=999  then  e{i)  =  .; 
if  ep(i)E=999  then  ep(i)  =  .; 
if  holder  e  1 .4  and  i  <54  then  do; 
e{i}  =  .; 
ep{i}  =  .; 
end; 

if  holder  =  9.7  and  i>53  and  moisture>13  then  do; 
e{i)  =  .: 
ep{i)  =  .: 
end; 
end; 
drop  i; 
run; 

proc  printto  new  print  =  'auQ93\aug93_5.prn'; 

options  ps  =  54; 
run; 

proc  print; 

var  date  material  holder  hidrnum  moisture  temp; 
run; 
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SAS  Code  for  Generating  Frequency  Plots 


•  program  "plotallf",  directory  "augSS” 

•  « 

$ 

« 

*  "powattn"  is  power  loss  -20log10(expM/skin  depth)}  in  dB/meter 

• 

*  "phasevel"  is  the  ph^e  velocity  »  radial  frequency/beta 

*  "sigma"  is  the  electrical  conductivity  in  mho/m  and  is  calculated 

•  by  the  formula:  epspp* radial  freq'free  space  permittivity 

• 

*  "losstan"  is  the  loss  tangent  defined  as  epspp/epsp 

libname  eps  'c:\sas_win\aug93\'; 
data  plt^tmp; 
set  eps.aug93_5; 
array  e{100}  el -el  00; 
array  ep{100}  epI-eplOO; 

array  f{100}  (.045,.078..1 1 1,.144..177..210,.243,.276..310,.343. 
.376,.409,.442,.475,.508..541,.574..607.,640..673,.706,.739,.773, 
.806,.839,.872,.905,.938,. 971, 1.004. 1.07, 1.1 03, 1.1 69. 1,235,1. 269, 

1 .335,1 .401 .1 .467,1 .533.1 .632,1 .698. 1 .765,1 .864,1 .963,2.029.2,1 28, 

2.228,2.36,2.459,2.558.2.691,2.823.2.955.3.087.3.253.3.418,3.550, 

3.749,3.914,4.112,4.31 1.4.509.4.708.4.939.5.171.5.435.5.667.5.964. 

6.229.6.526,6.857.7.188.7.519.7.882,8.246.8.643,9.073,9.503.9.933. 

10.429,10.925.11.454.1 1.983.12.545.13.173,13.802.14.43.15.124. 

15.852,16.612,17.406.18.233,19.126.20.019,20.978.22.003,23.061. 

24.152,25.310,26.5); 

if  materiar  =  'AIR'; 

drop  date  hidrnum; 

do  i  to  100; 

epsp  =  e{i); 

epspp=:ep(i}; 

if  epspp<0.2  then  epspp^O.; 
freq=:f{i}; 
powattn  =  .; 
phasevel  == .; 

sigma  =  epspp*2*3.141 59  "freq*  0.00885; 
losstan  s  epspp/epsp; 
if  epsp  =  999  then  losstan  =  .; 

phasevel  =  1  /(sqrt((epsp * (sqrt(  1  +  (epspp/epsp)  *  *  2)  +  1  ))/2)); 

if  epsp  =999  then  phasevel  = .; 

if  epspp=0  then  go  to  skp; 

if  epspp  =  .  or  epsp  =  .  then  go  to  skp; 

num  =  3e8/(freq  *  1  e9); 

powattn  =  num/(2*3.141 59*sqrt((epsp*(sqrt(1  +  (epspp/epsp)*  *2)-!  ))/2)); 
powattn  =-1  /powattn; 
powattn=-20*log10(exp(powattn)); 
skp:  tmp  =  temp; 
output; 
end; 

drop  el-el 00  epI-eplOO  fl-flOO  i  num; 

run; 

• . 

proc  sort; 

by  material  temp  descending  holder  moisture: 

run; 

•  • 

goptions  reset  =  all  device  =  hpljs3; 
titlel  h  =  2  '  '; 
titles  '  '  a=90  '  '  a  =  -90  ' 
footnote  h  =  6  ' 
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axisi  length- 13  cm  logbase  =  10  order  =  (.01  .1  1  10  100); 
axis2  length^  13  cm  logbasess  10  order  =  (.l  1  10  100); 
axis3  length  =  1 3  cm  logbase  =  10  order-(.001  .01  .1  1  10); 

axis4  iength  =  13  cm  order  =  (0  to  1  by  .1); 

axisS  length  =  13  cm  logbase  — 10  order  =  (.01  .1  1  10  100  1000  10000); 

symboll  1  =  1  i=jw  =  2; 

symbol2  1  =  2  i=j  w  =  2; 

symbols  1  =  3  i=j  w  =  2; 

symbolA  I  =  4  i  =  j  w  =  2; 

symbols  1=8  i=J  w  =  2;  , 

symbols  1  =  14  i=j  w  =  2; 

symbol?  I  =  33  i  =  j  w  =  2; 

symbols  I  =  35  i  =  j  w  =  2; 

symbols  1=41  i  =  j  w  =  2; 

symbollO  1=43  i=j  w  =  2; 

symboll  1  1  =  1  i=j  w  =  2; 

symboll  2  I  =  2  i=j  w  =  2; 

• . 

proc  gplot  data=plt_tmp; 
by  material  temp; 

plot  epsp*freq  =  moisture  /haxis  =  axisi  vaxis  =  axls2; 
label  epsp  =  'REAL  EPSILON'; 
label  freq  = 'FREQUENCY  ~  GHz'; 
label  temp  =  '  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE'; 
label  material  = 'SOIL'; 
run; 

». 

proc  gplot  data  =plt_tmp; 
by  material  temp; 

plot  epspp*freq  =  moisture  /haxis  =  axis1  vaxis  =  axis2; 
label  epspp  =  'IMAG  EPSILON'; 
label  freq  = 'FREQUENCY  ~  GHz'; 
label  temp  =  '  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE'; 
label  material  = 'SOIL'; 
run; 

•  • 

proc  gplot  data  =plt_tmp; 
by  material  temp; 

plot  sigma •freq  =  moisture  /haxis  =  axisi  vaxis  =  axisi; 
label  sigma  = 'CONDUCTIVITY  ~  mho/m'; 
label  freq  = 'FREQUENCY  ^  GHz'; 
label  temp  = '  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE'; 
label  material  = 'SOIL'; 
run; 

•  • 

proc  gplot  data  =plt_tmp; 
by  material  temp; 

plot  losstan*freq  =  moisture  /haxis  =  axisi  vaxis  =  axis3; 
label  losstan  =  'LOSS  TANGENT'; 
label  freq  = 'FREQUENCY  -  GHz'; 
label  temp  =  '  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE'; 
label  material  = 'SOIL'; 
run; 

«  • 

proc  gplot  data  =  plt_tmp; 
by  material  temp; 

plot  phasevel'freq  =  moisture  /haxis  =  axisi  vaxis  =  axis4; 
label  phaseve!  = 'NORMALIZED  PHASE  VELOCITY'; 
label  freq  = 'FREQUENCY  ~  GHz'; 
label  temp  =  '  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE'; 
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label  rnateria!  =s'S01L'; 

run; 

• . 

proc  gplot  data  =plt_tmp; 
by  material  temp; 

plot  powattn*freq  =  moisture  /haxis  =  axis1  vaxis  =  axis5; 
label  powattn  = 'POWER  ATTN  —  dB/m'; 
label  freq* 'FREQUENCY  -  GHz'; 
label  temp  =  '  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE'; 
label  material  =  'SOIL';  • 
run; 


I 
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SAS  Code  for  Generating  Moisture  Content  Plots 


•  program  "plotalfm",  directory  "augSS" 


*  "powattn"  is  power  loss  -201og10(expM/skin  depth))  in  dB/meter 

*  “phasover  is  the  phas^e  velocity  =  radial  frequency/beta 

* 

*  "sigma"  Is  the  electrical  conductivity  in  mho/m  and  is  calculated 

*  by  the  formula:  epspp*  radial  freq*frce  space  permittivity 

* 

•  "losstan"  is  the  loss  tangent  defined  as  epspp/epsp 

•  « 

libname  eps  'c:\sas_win\aug93\^' 
data  plt_tmp; 
sot  eps.aug93_5; 

array  f{100)  {.045,.078,.1 11  ,.144,.177,.210,.243,.276,.310,.343. 
.376,.409,.442,.475,.508,.541,.574,.607..640..673,.706„739„773, 
.806,.839,.872..905,.938,. 971 ,1 .004,1 .07,1 .103,1 .1 69,1 .235,1 .269, 

1 .335,1 .401 ,1 .467,1 .533,1 .632,1 .698,1 .765,1 .864,1 .963,2.029,2.1 28, 
2.228,2.36,2.459,2.558,2.691,2.823,2.955,3.087,3.253,3.418,3.550, 
3.749,3.914,4.112,4.31 1,4.509,4.708,4.939,5.171,5.435,5.667,5.964, 
6.229,6.526,6.857,7.188,7.519,7.882,8.246.8.643,9.073,9.503,9.933, 
10.429,10.925,11.454,1 1.983,12.545,13.173,13.802,14.43,15.124. 
15.852,16.612,17.406,18.233,19.126.20.019,20.978,22.003.23.061. 
24.152,25.310,26.5); 
if  materiar  =  'AIR'; 
drop  date  holder  hldrnum  temp; 
epsp  =  e2; 
epspp  =  ep2; 
f  req  =  0.078; 
powattn  3= .; 
phasevels.; 

sigma  =  epspp  *2  *3.1 41 59 'freq  *0.00885; 
losstan  =  epspp/epsp; 
if  epsp  =  999  then  losstan-.; 

phasevel  =  1  /(sqrtdepsp *  (sqrt(l  +  (epspp/epsp) *  * 2)  + 1  ))/2)); 

if  epsp  =  999  then  phasevel  =  .; 

if  epspp  =0  then  go  to  skp; 

if  epspp-.  or  epsp  =  .  then  go  to  skp; 

num  =  3e8/(f  req  "le9); 

powattn  =  num/(2*  3.1 41 59  *sqrt{(epsp*(sqrt{1  +  (epspp/epsp)*  *2)- 1  ))/2)); 
po  wattn  =  - 1  /po  wattn; 
powattn  = -20  *log10(exp(powattn)); 
skp:  tmp=temp; 
output; 

drop  el -el  00  epI-eplOO  fl-flOO  I  num; 
run; 

«. 

proc  sort; 

by  tmp  material  moisture; 
run; 

•  • 

goptions  reset  =  all  device  =  hpljs3; 
titlel  h  =  2 

title2  j  =  c  h  =  2  f  =  swiss  '78  MHz'; 
title3  '  '  a  =  90  '  '  a  =  -90  '  '; 
footnote  h  =  6  ' 

axisi  length  =  13  cm  order  =  (0  to  100  by  10); 

axls2  length  =  13  cm  order  =  (0  to  100  by  10); 

axis3  length  =  1 3  cm  logbase  =  10  order  =  (.001  .01  .1  1  10); 

axis4  length  =  13  cm  order  =  (0  to  1  by  .1  ); 
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axis5  length  s  1 3  cm  logbase  « 10  order  =  (.01  .1  1  10  100  1000  10000); 

axis6  length  =  1 3  cm  logbase  s*  1 0  order  =  (.01  .1  1  10  100); 

symbol  1  v  =  A; 

symbol2  v  =  B; 

symbols  v  =  C; 

8ymbo)4  v  =  D; 
symbols  v  =  E; 
symbols  v  =  F; 
symbol?  v  =  G; 
symbols  v  =  H;  , 

symbols  v=l;  • 

symbollO  v=J; 
symbol  11  v  =  K; 
symbon2  v=L; 
symbollS  v=W; 

« . 

proc  gplot  data  =  plt_tmp; 
by  tmp; 

plot  epsp*moistures=material  /haxis  =  axis1  vaxissaxis2; 

.  label  epsp  =  'REAL  EPSILON'; 
label  freq  = 'FREQUENCY  -  GHz'; 
label  temps'  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE  (VOLUMETRIC)'; 
label  material  = 'SOIL'; 
run; 

• . 

proc  gplot  data  =plt_tmp; 
by  tmp; 

plot  epspp •  moisture  =  material  /haxis  =  axis1  vaxiss8xis2; 
label  epspps'lMAG  EPSILON'; 
label  freq  = 'FREQUENCY  GHz'; 
label  temp  = '  TEMPERATURE'; 
label  moisture  *'%  MOISTURE  (VOLUMETRIC)'; 
label  material  =  'SOIL'; 
run; 

«• 

i 

proc  gplot  data  =  plt_tmp; 
by  tmp; 

plot  sigma  •moistures  material  /haxis  =  axisl  vaxis  =  axis6; 
label  sigma  = 'CONDUCTIVITY  ~  mho/m'; 
label  freq  = 'FREQUENCY  -  GHz'; 
label  temp  =  '  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE  (VOLUMETRIC)'; 
label  material  'SOIL'; 
run; 

• . 

proc  gplot  data  =  plt_tmp; 
by  tmp; 

plot  losstan •moisture  ~ material  /haxis  s^axisl  vaxis^^axisS; 
label  losstan  =  'LOSS  TANGENT'; 
label  freq  = 'FREQUENCY  ~  GHz'; 
label  temp  s'  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE  (VOLUMETRIC)'; 
label  material  =  'SOIL'; 
run; 

». 

proc  gplot  data  =  plt__tmp; 
by  tmp; 

plot  phasevel* moisture  =  material  /haxis  =  axis1  vaxis :saxis4; 
label  phasevel « 'NORMALIZED  PHASE  VELOCITY'; 
label  freq  = 'FREQUENCY  ~  GHz'; 
label  temp  =  '  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE  (VOLUMETRIC)'; 
label  material  = 'SOIL'; 
run; 
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proc  opiot  data  =  plt__tmp; 
by  tmp; 

plot  po wattn* moisture  s=  material  /haxis  =  axis1  vaxis  =  axis5; 
label  powattn  = 'POWER  ATTN  ~  dB/m'; 
label  freq  = 'FREQUENCY  —  GHz'; 
label  temp  = '  TEMPERATURE'; 
label  moisture  =  '%  MOISTURE  (VOLUMETRIC)'; 
label  material  'SOIL'; 
run; 
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